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The Mechanics of Non-Linear Materials with Memory 


Part LEE 


A. E. GREEN & R. S. RIVLIN 


1. Introduction 


The mechanics of non-linear materials with memory has been discussed in 
previous papers [/, 2] and an earlier report [3]*. In these papers, it was assumed 
that, in a fixed coordinate system, the stress at a point of the material at any 
instant of time is determined by the deformation gradients at that instant and 
at previous instants. The limitations imposed on the form of the constitutive 
equation by the consideration that it is unaltered by a simultaneous rotation 
of the body and reference system are determined. In the present paper, we 
assume initially a constitutive equation in the form of implicit relations between 
the stress, its time derivatives and gradients of displacement, velocity, accelera- 
tion, etc. at a number of instants of time in some time interval. The limitations 
on the form of the constitutive equation resulting from the fact that it is un- 
altered by the imposition on the body of an additional arbitrary angular velocity, 
acceleration, etc. are discussed. 

In §§ 6 and 7, we make the assumption that the stress in a fixed coordinate 
system at any instant of time depends explicitly on the gradients of the displace- 
ment, velocity, acceleration, efc. at that instant and previous instants. The 
limitations imposed on the form of such a constitutive equation by the considera- 
tion that it is unaltered by a simultaneous rotation of the body and the reference 
system are determined. 

The remainder of the paper is concerned with the relation between such a 
constitutive equation and constitutive equations based. on apparently more 
limited initial assumptions. It is shown in §6 that if the dependence of the 
stress at any instant on the deformation history at preceding instants is suffi- 
ciently smooth, and if the motion starts from rest sufficiently smoothly, then 
there is no substantial loss of generality in assuming that the stress depends on 
the displacement gradients only at preceding instants and on the gradients of 
displacement, velocity, acceleration, etc. at the instant of measurement of the 
stress. In §8, it is shown how various limitations in our initial assumption 
regarding the arguments on which the stress depends is reflected in limitations 
on the form of our resulting constitutive equation. 

In §9 it is shown how particular dependence of the stress on the gradients 
of displacement, velocity, acceleration efc. may be formally included in functional 
dependence of the stress on the displacement gradients only by the use of Dirac 


* The subject has also been discussed by Nort (4, 4]. 
Arch. Rational Mech. Anal., Vol. 4 26 
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delta functions and their derivatives. Finally, in §10, the relation of memory 
theories of the type discussed in the present paper to a previous theory of RIVLIN 
& ERICKSEN [6] is discussed. 


2. Kinematic Preliminaries 


Let every particle of a three-dimensional body be at rest, prior to time tT=0, 
relative to a fixed rectangular Cartesian coordinate system x. Let us suppose 
that after zero time the body is deformed continuously so that a typical particle 
P, of the body moves to P at time ¢ and that at any time t (0StS?) the co- 
ordinates of the particle in the system are x,;(t). We shall employ the notation 


X;=%,(0), j= %,(8)- (2.1) 


We consider that the functions x;(t) are specified as continuous functions of X; 
and t with continuous partial spatial derivatives up to any required order, 
except possibly at singular points, lines and surfaces. We then have 


x(t) =%,(X;,7), %;,=%;(X;,2). (2.2) 
If this deformation is to be possible in a real material then 
Ox; (t) 
“ax, | 77 9 (2.3) 


The components of velocity of the point P in our fixed Cartesian frame are 
denoted by v! (rt) =v,(t) and we use the notation v,(f)=v;. Thus 


Dx; (7) 


vy (rt) =;(t) = (2.4) 


where D/D rt denotes differentiation with respect to t holding convected or material 
coordinates X; fixed. The components of acceleration v{*)(z) are then given by 


Dv (x) ev (t) 
H(t) = = ot (2) 


1 =D, 


avy? (t) 
OXm(t) 7 


(2.5) 
where 0/07 denotes differentiation with respect to time keeping x,(t) constant. 


More generally, the (n —1)" acceleration components may be expressed as 


Dof—» (x) ) 

(n) =. v (1) é (n—1) 

ON Ga 2 4 | Cat 

(7) Dt ae tm (7) OX, (T) 2 (z), 


vi w(t), vl (cz) = x, (1). 


(2.6) 


We now consider another deformation of the body which differs from that 
described by (2.2) only by a superposed arbitrary rigid body rotation. Then the 


coordinates %,;(t) of the particle P, at time t in the coordinate system 4%, are 
given by = 
X;(t) =a,,(t) x,(t), (2.7) 


where a, ;(t) are any set of continuous functions of t defining a rigid body rotation 
and are therefore subject to the conditions 


4; ,(T) a;,{T) =, A(T) a, ;(T) Pp Ons | a;;(z)| =a (2.8) 
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We denote the components of velocity, acceleration, ..., (n —1)' acceleration, 
in our fixed frame, due to the resultant of the deformations (2.2) and (2.7) by 
vf” (t) =9;(r), Df (zt), ..., J (rt) and use the notation 9%) (t) =o") (s=4, 2,..., n). 
Then 
Dy) (t) =, 8f"(t) = %,(z), (2.9) 
ROIS — le) ere 
From (2.6), (2.7) and (2.9) we have 


aa IDEs Da;; 
7; (c) = PEO) — a,,(c) vj(c) + 8A a, (0) 
Daim (2) 7a alt 
= a;;(T) 4 ( Var IB jm (T) X; (7) 
Also, from (2.8), 
D dim (T) Dajm(t) 
3 ar A; (T) a Gim (t) ep pow a7 0, (2,44) 
so that «) (zt), defined by 
IDiGigep D Aj y(t 
aft) (x) fae A; m(T) oo ao oe A; m(T) aoa ? (2.12) 
is skew symmetric and we may define a vector Q")(t) by the equations 
aly? (7) = Cin Qh (t),- 20)? (t) = Cn ig OI (t), (2.13) 


where e;;,=1 or —1, accordingly as 7,7, is an even or odd permutation of 
1, 2,3 and e;,,=0 otherwise. With the help of (2.12) and (2.13) equation (2.10) 
becomes ss ~ 

U,(t) = 4,,(t) v,(t) + €:n¢ 2h) (7) X;(T) (2.44) 
and therefore Q(z) is the component along the x; axis of the angular velocity 
of the superposed rigid body motion at time t. The angular acceleration, second 
angular acceleration, ... of the rigid body rotation at time t are 


Ds-2 QW 
Oe ee a 
IDS 
1 (s) (2.15) 
== 5 C354 %h; (T)s 
where = 
F DEV oh) 
ab) (t) = (2.16) 
If we substitute (2.12) into (2.16) we obtain 
Sa 9) 
(s) a a Ajm(t) D> 4a; (7) 
a; (t) ou q Dra Dt—4 
and hence 
Seay q q 
IDE Giga (2) pS) B. Si “4 D Aj m(T) D Qj m (T) 
A; »,(T) Dt aa Cari (z) ya | q Dt Dea 
% q=1 Ate one (2.17) 
Mey Meena NTL Ses D1 ajm(t ain (7 
= of) e) — 3 (9) (tom) \ ant) 
By successive application of (2.17) we see that 
LDSdim 
aye ae (2.18) 


26* 
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may be expressed as the sum of «)(r) anda polynomial in «{(z), #{?(z),..-, an; (T) ; 
and we recall from (2.15) that the (s—1)' angular acceleration of the rigid 
body motion at time 7 can be expressed in terms of the three non-zero components 
of o{$) (1). 
3. Tensors Dependent on Rotation 
From (2.10) we have 


a5, (t) _ Bure) (ay Daim) 
0%; (t) <— ir (T) is | rae oe jm (T) Dt G 1) 
Ov, (T) D ajm(t) 
= 4;,(T) el ) Bx, (4) ik inl ) , De ) 
so that 
Daj,;(t) 80; (t) OVm (7) 
De ei) oe ep ae El eta 6.2) 
“eek 80 », (T) 6; (T) 
my a ) 0%; (7) + jm (T) on, 1@) 
Using the notation 
x 4 (00;(T) 60; (7) ) 
wij(t) =O (0) = 3 Lata) — (3.3) 
o,,(t) =o(=- IND ay See . 
bee a Gags 2 ae 0%;(t) I’ 
we see, with the help of (3.1) and (2.12), that 
@;;(T) = a;,(t) a;,(T) W,s (z) a ox) (T) z (3.4) 


If, at a particular time t=1,, the body is rotating but occupies the same in- 
stantaneous position at this time as that given by (2.2), we have a;;(t,) =6;,; and 
ij (Tx) = Oi; (Tx) + OY (Ta) - (3.5) 


Hence @,;,(t,) is not independent of the angular velocity of a superposed rigid 
body rotation at time f,,.. 


More generally suppose 


(s) ae du,” (rz) dv} (z) 
7; (7) 2 { 0x; (T) x; (7) \, 3 6) 
a (2) = 1 [2 ) _ 2) 
ti 2g { Om;(%) "0%; (z) } 
for s=1, 2,.... From (2.6) to (2.9) we have 
a(s 1S Wea m 
BAP ere seg 
00,” (z) = [s D*Taim(t) Ovi) (z) ie 
S(t 
Oa (ec) NG ie ID ig Oxy (Tt) ? 
where v!(t)=x,(z), 0 (t) =X, (z), and therefore 
BE} 2) = amt) nt) f(t) +5 {apn(a) BAO) — ay (0) PSD 4 
a a 8) 
Mens D® 4 aim (et) BvP (7) D*~4 aj (7) dum (tT) ) 2 
. 2 Dig) (eal Dt? Ax, (z) 4:4(T) De Oxp(t) J 
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If, at a particular time t=T1,,, the body is rotating but occupies the same instan- 
taneous position at this time as that given by (2.2), we have a;,(t,)=06,; and 


alt DS Sep 
o (t) = of) (T,) +{ ij (7a) 2 Hie} | 


ees a 3.9) 
as i ee ee ROU ta) DE a yh Cea) OO (ee al} | 
NG. 2 Bee OX; (T;) ID 0x; (tT) J - 
q=1 7 a 


Also, from (2.47) and (2.18) with 4 j(Ta) = Oe» the quantities D*a;;(t,)/D tT; 
(s=1, 2,...) can be expressed as the sum of «{*)(z,) and a polynomial in «!?) (t,) 
ae Oe (t,). Hence @)(t,) is not equal to of (x ,) but depends also on the 
three components of the (s—1)™ angular acceleration of the superposed rigid 


body motion at time t,, and on lower order angular accelerations. 


tA 


4. Tensors Independent of Rotation 


We denote the components of stress in the coordinate system x; at the particle 
F, at time t by o;,;(t) and employ the notation 


0;;—6;;(t). (4.1) 


Now consider a deformation given by (2.2) and an arbitrary superposed rigid 
body rotation (2.7). The stress components at time 7 in the coordinate system ~; 
due to the resultant of the deformations described by (2.2) and (2.7) are denoted 
by o;,;(t). The stress components at the particle J) at time t, referred to a moving 
rectangular Cartesian coordinate system which coincides with the fixed system 
x, at time t=O and rotates with the rigid body motion described by (2.7), are 
still o;;(t). Hence, in view of the tensor character of stress components 


6, ;(T) a, a;,(T) a; ,(T) 0,5(T) D (4.2) 
From (4.2) we have 
D&E tj D rs D a § D ir 
a A a fe) fo) a oh. 


With the help of (3.2) this equation yields 


Dar (Bink 0;(t) = 00; (T) 
aie = Op (t) OF» (7) Om; (T) Xm (7) (4 3) 
1, (T dv, j 
— 1G) js(T) ee Sea ) 7G ee ec (t) a rl 
and 
D6;;(t) = 0Vm (T) = 20m (T) 
De + Pm() azcay tOmi\) a (4.4) 
= aip(t) aj5(2) [PF + 0, (1) at) + On s(t) mr |- 
Also, by adding (4.3) and (4.4), and using (3.3), we have 
IOGiRy = 
ST ENO T) + Gay (t T) Omi (T) (4.5) 


Doy5(T T) 
Dt 


ASG aE ) Wyn 5 (T) Gat) Wy (T) e 
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If we now introduce the notation 


D oj; (tT) OUm (7) OVy (T) 
of (1) = ae ot Gym) “8x; (1) + Gq; (T) Qx,(t) ’ 
D0;;(t) du; (Tt) 6u;(T) 
or (time ae ment ae AG AT) Diy (t) ” (4.6) 


The results (4.6) and (4.7) have been derived directly from (4.2) using also 
(3.3) and are valid whether or not o;,(t) is the stress tensor, provided it satisfies 
an equation of the form (4.2). Hence by repeated differentiation using (4.7), 
we find that 


6) (1) = Aj (T) An (T) Own (T) 
Bip CT) Ow (O) en On Ce) (4.8) 
QO Ga) DSC pe) 
a A; ,(T) A;(T) aa 5) 
for s=1, 2,..., where 
Do§) (z) OUm (T) Um (tT) 
(s+1) Me NID (s) He (s)_ ee 
07; (Z) Dt ali GAAe) x; (T) a Omi (T) 0x; (T) 
Do; (t) uv; (t) 6; (tT) 
* (s+1) ee * (s) 7 ey) $ ‘ 
07; ( ) Dt Vm (t) OXm(T) m7 ( ) O%m (7) ) (4 9) 
DoT a(t) D {¢ Bo;;(t) DB 65m(t) DB Om; (T) 
Rett ~ \ Al De } oe Omi) Vga Omit) 


and 


We now obtain an alternative expression for o{)(t) which will be of use 
later. From (2.7), (2.8) and (4.2) we see that 


G,j(t) 4, (1) d,(t) =0;;(t) dx,(t) dx;(x) 
so that f(r), defined by 

B(t) =0,,(t) dx,(t) dx,(t) =t,,,(t) dX, aX,,, (4.10) 
is independent of rigid body rotations superposed on the original motion, where 


tnt), defined by ie ne 
1 y ‘s 


are functions independent of rigid body rotations, for all values of t. This is 
also evident from the fact that ¢,,,,(t) are the covariant components of the stress 
tensor referred to curvilinear coordinates X,, in the deformed body at time T. 


Non-Linear Materials with Memory. III 393 


From (4.10) we have 


DB(t) — _Do;;(z) Didx; Didx; 
De pe alt) (2) $05j(0) SE ax, (2) to,;(0) dx, (2) 
_ Dtnn() 
=n dX, aX,,. 
Since 
D dv; (7) 
Dr dx; (tT) — dv, (t) —— 0x; (7) i(7) 
we have, using (4.6), 
D D 
a Ope dae aX, dX, (4.12) 
Repeated differentiation of (4.12) gives 
IDP DSt n 
ee ee oy a (4.13) 
Returning to (4.10) we see that 
DSGRa) a TsO ost) Sr gy Dida (2s) Di de;(zy 
Dis ic ()-se5 ys (;) Dita’ Dat ; eee 
Since 
D' dx; (zt) dv (7) 
2g 1 Sy) (k) aa 4 
Dx (7) 0x4 (T) Bole) 


and dv (t)/dx,;(t) =6;;, equations (4.14) and (4.13) yield the alternative formula 


6) Do omnl®) 149) 285, (x) ouh (2) 
(3) Digs > (i) et 0x; (7) 


0 
_ Droute) i D4 6mn(t) < ang (x) auth (x)) 1? 
— Dt cp lel i ae 23 0x; (T) a 
for o) “llce ). Also, from (4.13), 
DS tin n (T) Bees) 0x, (t) Ox; (7) 
Diss 2) \t) earls (4.16) 


and these functions are completely independent of any superposed rigid body 
rotations. 


Again, if 


ci a SGM > Coa) OS 
UNV SO Nt) aren az, ey Cantey Oa) ae 
then ?””(t) are independent of rigid body rotations for all time t and we may 
show that 


DOE ie. * (s) OXm OXy 
Saree 0 ia ) Bale) Ox; (tx) aA) 
The functions 7?””(z) are contravariant components of the stress tensor at time tT 
referred to curvilinear coordinates X,, in the deformed body at this time. 
When o;;(t), instead of denoting the stress tensor, has the special value 6,; 
we denote an corresponding values of o{)(t) and ¢,,,,(t) by A) (t) and gm,(7) 
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respectively. Hence, from (4.6), (4.8), (4-9), (4.15) and (4.16), we have 


ae Ale) = \ 
2A,,(t) oi A ] (T) x; (T) ax; (tT) ) v7 


A 
OX; (T) OG o)ee (4.19) 
v(t) dv (zt), (s\ avs (tz) av (zx) 

“Oaj(t) saa + DL @x;(t)  Ox;(t) ’ 
as L(t) =aim(t) an (t) AM4(2), 
°Smn(t) 2 gery 22), Oa) 
MT ee AN (*) 3x aX, (4.20) 
OX; OX; 
Ent) = ee 


When o;,(t) =0,;; we observe from (4.6) that 
Di; 
Bt 
This relation is of importance when we express our results in terms of a general 
fixed curvilinear system of coordinates. 
If at time t=T, the body has a superposed rigid body rotation but occupies 
the same instantaneous position as that given by (4.2) (with t=T,) we see from 
(4.8) and (4.19) that 


Gey (Ea) 208) (eg) ete (alae A ee a eee dee ee (4.22) 


=; (4.21) 


5. General Constitutive Equations 
We consider general constitutive relations of the form* 


j Ore (Tal) Oeste) “Con ce) dup” (Tx) 
Wh (OX Gg OA Som SaOaa\ Gala. eee) 5.) 
D 6p, (Tz) DOE O sal Gen) 
pq\°a pq\ta 
oO ie —————,... , | = 
Pye a) Dtz 5) , Dim O; 


where t, («=0, 1, 2,..., N) denote N+ distinct instants of time between Tt=0 
and t=?, and where t)=?. Also we assume /;;=/;;. From (4.15) we see that 
D"6y4(t,)/D ty can be expressed as the sum of o§”)(r,) and a polynomial in 


DO Op {TaD te +22; DO pq tall D Ten God alee Oe ee/ONgltx)s5. one, (Ug erete te 


with similar results obtained by replacing m by m— 1, m— 2, ..., Hence, equations 
(5.4) can be written in the different form 
a OXp (Tq)  OXp(t,) “Ovs?(z,) dv}? (Ty) 
eo Pe BCG. SNP AR (5.2) 
Opolta)s Oe Coens of) (r4)| =0 (n=m), 
pl Saeed, Bante 20 Era) Bere 
v7) OX, ? Ox, 5) 8x4 (T) 1G) Ou Oay eae” (s 3) 
Opq (Ta), Thy (Ta) ---, OM (t,)|=0 (nm). 


cof fro) order to specify the constitutive equation completely, we must in general 
state appropriate initial conditions on the differential equations (6x1); 
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Moreover if (5.1) is a polynomial in its arguments so also are (5.2) and (5.3). 
Without loss of generality we may consider equations of the form (5.2) where 
n,m are integers such that n=m. 

With the help of (3.6) and (4.19) we may express 2év{" (r aioe (q(T) as the 
sum of 208) (r J+A pia, eamd a eae. in the gradients vy) (T.)/OXq (Tq) for 
f—7—_A, n= 2, ae Similarly 2005" ”(r,)/0x,(t,) can be replaced by fe sum 
Choe. (r 44g Y(t,) and a Bolynouiial in the gradients dv) (r,)/Ox,(t_) for 
y=n—2, n—}3,...,1, and so on. Hence, by successive replacement, we can 
reduce (5.2) to fe different form 


O% p(T.)  O%p(T.) 
fii] gx aq? Abela) ABN (ta); (5.4) 
Opal Cea salem Oey MEVA@ Se) (Ta). Oe (ey) (== Ou (0 =m), 


If (5.2) is a polynomial form we see that (5.4) is also a polynomial. Hence, if 
(5.1) is a polynomial form so also is (5.4). On the other hand if (5.1) is a continuous 
function of its arguments with continuous derivatives up to some order, so also 
s (5.4). 

We now consider a motion of the body which differs from that described 
by (2.2) by a superposed rigid rotation. Such a motion is given by (2.7) and 
(2.2). Now suppose that at the times t, («=0, 1,..., N) the body is in the same 
instantaneous position as that given by (2.2) at these times, but that it has 
superposed angular velocity, angular acceleration, ..., (7 —1)™ angular accele- 
ration on the motion of every point of the body. Hence a; ;(t,) =6;;(«=0, 1,...,N) 
and in view of (4.22) equations (5.4) are replaced by 


Cee) OX (T.) 
Ai pase (1) () (7 )- 
hij aX, ? Ox, , ree (ee Were ral a3 ae 
Opq(Ta) » O57 (Fa), hes Oe CANT ec cope (t, )|=0 (n=m), 
where @$}(Tx),---,@pq(T) are given by (3.9). If the form of the constitutive 


equations is to be independent of the three components of the (7 — 1)" angular 
acceleration of the superposed rigid body rotation, at each of the times t,, then 
since the (7 —1)'" angular accelerations are independent at each of these times 
wy") (t,) («#=0,1,...,.N) cannot occur explicitly in (5.4). If wl") (z,) is removed 
from (5.4) and @{?(zr,) is removed from (5.5) then the form of the constitutive 
equations will be independent of the (m— 2)" angular acceleration of the super- 
ae rigid body rotation at each of the times Tt, if Con Y(r,) does not occur 

n (5.4). Hence, by successive reductions of this type, we see that the equations 
ee will be independent of any angular velocity, angular acceleration, ..., (7 — A\ 2 
angular acceleration of the rigid body rotation which may be superposed at 
times t,, the body occupying the same instantaneous positions at these times, 
if the equations (5.4) reduce to 


edu) Chale - i ps 
fii: i M a A (14), AY (tg) pg (Ta); 09) (tq), «++ 0%) (4) | = 0, 
(n=m). (5.6) 


So far our considerations have been purely formal. We shall now show how 
they can be applied in particular cases. 
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6. Functional Dependence of Stress on Kinematic Variables 
Let us now suppose that the stress components o;,; at time ¢ are functionals 
of the gradients 8x,(t)/@X,, Ox,(r)/6x,, Ovp?(t)/Ox,(z), dup (a) (ae \ulsena 
du” (r)/Ox,(r), thus: ? 
(6.1) 


We first replace the functional F;,; by a polynomial function /;; of the values 
of the arguments at N-+1 discrete instants of time t, (=O, 1, 2,..., N), where 
T=t, and obtain 

Oty (Ga): Op (tq) Op Fa) avg” (r.) 


, d ’ aa) a (6.2) 
0X4 OXg OXq (T.) 0% (Tx) 


O4=Hi; 


Following the argument of § 5, we see that the relation (6.2) must be expressible 
in the form 


0% p (Tx) OXp (Tx) (1) (n) 
Sees ALY), ASY(D)), (6.3) 


O;=hi; 


where the polynomial function /;; is not in general the same as that in (6.2). 
We now let N-+1-— oo, so that the function /;; in (6.3) becomes a functional. 
We then have . 


at O%p(t) O%p(T) 4a (n) 
04 =F] es pe ABO), «2s A (|, (6.4) 
7—0 


where the functional F;, is not in general the same as that in (6.1) *. From (4.20), 
we have 
DS ginn (7) on 0x; (T) 0x; (7) 
S a v7 
and i Wikia a sci (6.5) 
m n Bmn(T) 
0x,;(t) Ox;(t) Dice i 


We note also that 


OXg a ax, OX, (6.6) 
and 

Ox, as pit (igs atee9) 

Big. ah Ree ae (6.7) 


fo 


From (6.6) and (6.7) we see that any constitutive equation of the form (6.4) is 
equivalent to a constitutive equation of the form 


t 


ax» (1) : 
oi =F aX, ARID), AGHA). (6.8) 


From (6.5), we note that (6.8) may be expressed in the form 


(6.9) 


OX, er 
™=0 


78 INKS shall use £;; to denote a functional generically in many places below. They 
will not, in general, all denote the same functionals of the indicated arguments. 
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We note that if in (6.4) F,, is a continuous functional* of the argument functions 
in the interval O<t<t and particularly dependent on their values at t=1, 
then so, in general, is the functional F,, in (6.9). We shall assume that the 
particular dependence on the values of the argument functions at t=¢ is suf- 
ficiently regular so that F,, in (6.9) may be expressed as a polynomial in 0% (0X4, 
Dg,,/Dt,..., D”g,,/Dt", the coefficients in this polynomial being continuous 
functionals of the argument functions over the range 0<1<t, and the notation 
Spqg=8pq(t) is used. By a procedure similar to that used in a slightly different 
context in previous papers [J, 2, 3], each of the coefficient functionals may be 
uniformly approximated by the sum of a number of terms of the form 


t ¢ t 


| 6p,(ts) Op, (ta) 
| i x “| Leen ile nace faaal T1, Ve, -5-, Tr) ax, ae 2 
06 66 : a 
D 8pn9% (TK) er D? pz qr (TL) ri D" gman (TR) ab be (6.10) 
Dtx Dt Dt, PLY Om? Ue 


where the y’s are continuous functions of their arguments. 


Now 


t 
Ds gp "gd (Ty) 
f etteivupaan (ls Th Varese te) — EO dy 
tN 
0 


t 

D Ds" Spy qv (tN) 
=. {sim beanl Ty, Ta, +++, TR) a) et = \ dry 
0 


t 


DY ai eR Ga as Ph n>, © py an \EN) 
dies ww\TN) g ; 
/ Dty Dry > Bae 
0 
IDS (t ) t 
ea ‘ee Paty l N 
Heese eo baer! 2 io TR) Diy? ty =0 


t 
D Hiji,j,.pran(h Tas Tas +--+ FR) Ds gon ay (tn) 
Dtyx Dty} 


dty. 
0 


We assume that the motion starts sufficiently smoothly so that when tT=0*, 
vs) (rz) (s=1,2,...,2—1) vanish and hence D*g,,(t)/Dr° (s=1,2,...,%—1) 
vanish. Then, provided the functions involved are sufficiently regular so that 
the differentiations may be carried out, we see that we may, without further 
loss of generality, replace the functionals (6.10) by functionals which are con- 
tinuous functionals of 8x,(t)/@X, and polynomial functions of 6x,/0X,, Dg, ,/Dt, 

.., D"g, Dt" and hence of 0x,/9X,, A$q,-..,Ap,- Thus, in (6.8), we may 
take o;,; to be a polynomial function of 0% (OX g, ise ee A” and a continuous 
functional of 0x,(t)/0X,, so that 

t 
a ee eee) Awl 


Om) v7 ax, > 6X, , PQ See 
t=0 


(6.12) 


* Throughout this paper, where we describe a functional as a continuous func- 
tional, we shall mean that it is a continuous functional over the compact aggregate 
of the argument functions which are continuous in the range OStSt. 
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Similarly, in (6.9), we may take o,, to be a polynomial function of 0x,/0X,, 
Dg, (Dt, ..., D®g,,{D and a Saunn tana functional of 6x,(t)/0X,, so that 


ax @ lex ODE Dg 
La= eat p ee pee oLde en 6.1 
055 =F; TORRE GC, 2 SDPSY 4 Tae (943) 
=. 


In the next section we shall consider constitutive equations of the types (6.12) 
and (6.13) further. 

The results (6.12) and (6.13) may be obtained by a slightly different procedure. 
As before, we see that we may without loss of generality omit the argument 
functions 6x,(t)/0x, in (6.1). We then assume that in (6.1) o;; is a continuous 
functional of the argument functions in the range 0<1t<¢ and is particularly 
dependent on their values at t=?. Using the relation (6.7) it eee that o;; 
is a continuous functional of 0x,(t)/dX,, du, (OX, =: rou Hes )/oX, in the 
range 0<t<¢?# and an ordinary function of the values és these functions at 
t=t. We assume sufficient regularity so that o;; may be uniformly approximated 
by a polynomial in 0% (OX g, G0, 10K, ee aul aX,, the coefficients In which 
are continuous functionals of 0x,(t)/0X,, dv) (t)/OX,, ...; dvy (x)/OX,. A typical 
term in the expression for o;,; 1s then 


gu) aul) a 
aX; OX; a? fe ee ARIR PiU psas (ts T, Tz, +++) Ts) X 
6.14 
aux) Buy?! (zk COUs as) via 
Xx “Ps Be oP Ce I ag fe 
OX, 0X4, 0X5 a 
where the kernels are continuous functions of their Dee ee Oy Hay 2. Ole 
and f,, By,..-, By are chosen from the integers 0, 1, 2,..., m and v" denotes x;. 
We note that 
ay (by) (ty) 
b N 
[xis Peco aglts 73-74 Ts) a dt 
(By —1) (6.4 5) 


(Bv—1) (4. : 
Wane epee OUpy (tw) | i D Kijigin P88 OV by (tw) dtn 
LJ Lee 8 r,) axe = ny : D Ty OXoy N 
By repeated use of relations of this type it is seen that o; ; may be get 
approximated by a polynomial function of 0x,/0X,, dv} aX, 1 Ooe 8X, the 
coefficients in which are continuous functionals of OX» ® \/aX,. rayne 
t 

Ae Tt) 0%» du? dup” 

H de OX petra deny OX, eae One, 


t™=0 


(6.16) 


Using the relation (6.7) it is seen that o;; May be PO HED, approximated by a 
polynomial function of 0x,/0X,, dul HCE - LD 1 Oka, the coefficients in which 
are continuous functionals of 0x, (t)/OX,. isn, 


ies Yay (tae aae ee ebance 6 
vq Oi) OX, 5} aX, ) ax, Shy UROL) Ox, . ( 17) 
t™=0 


In a manner similar to that used in deriving (6.12) from (6.1) we can show that 
(6.17) must be expressible in the form (6.12). 
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7. Further Restrictions on the Constitutive Equation 


We take as our starting point a constitutive equation of the form (6.13) 
and investigate the further restrictions which may be imposed on it, short of 
assuming that the material has some symmetry. 


From (4.17), we have with the notation ¢”” (1) =/”"” 


ne, OX; Ox; mn 
NEES Geo els 
and (7.1) 
iy ages yas wae 
i Forte xe: 
It follows from (6.13), (7.1) and (6.7) that 
t 
‘j Oxp(t)  O%, Dg D” g 
Ope ge p eee Ps Spq pq 
f ey OX e Oak Dt pa 2H ott” (7.2) 
™=0 


where £;,; denotes a different dependence on the arguments from that in (6.13). 
We replace the functional dependence of #’’ on 0x,(t)/OX, by polynomial 
dependence /;; on the values of the functions at N-+1 discrete instants of time 
GAle—=O, 15/2, %60., N+ 4)) with ot; thus: 
; Oxp(t,) De D"g 
er b\To 5pq Pg 
t fail xe oD SEDI (7.3) 


We next consider a deformation of the body which differs from that described 
by (2.2) only by an arbitrary rigid-body rotation (2.7) in a small time interval 
about t,. From § 4, we see that ios DgpDi, ..., D’ gp {Dt and 0%, (t,)/eX,, 
(t,=-T,), are unaltered by this rotation, so that using (7.2) we have 


i O% p(t) O%p(Tx) D8pq ae 
lah LS al uals ee (7.4) 
sa X(t) Xp (Tx) D8pq =e] 
ij aX, 5 aX, , Dt pottery Di” 


Thus, f;; is a polynomial scalar invariant under proper orthogonal transformations 
of the vectors 0%,(T,)/OX1, 9%» (T,)/OX.» and Ox» (t,)/OX3. Hence we see that hij 
must be expressible as a polynomial in g,,(t,) and [g(t,)]*. Repeating this 
procedure for each value of t, in turn, we see that 


47 1 Dg, IDierg 
2 un eee dala ee see eee (7.5) 
Now, allowing N +1 o9, we obtain 
t 
ij Dg DEE, 
Gi = il tpal@))- Ceo pees o> pepe (7.6) 


t=0 


From (7.1) and (7.6), we obtain 
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Using (6.5), we may re-write (7.7) in the form 


t 
Ox, Ox; ; , Oxy OX (1) OX, OX5 A], (7.8) 
[i= Ox, OX, 1” Fane POH gars gmt e ts 1 ORR eH 2 50 
se 


In (7.6), (7.7) and (7.8), F, is a continuous functional of the argument functions 
q(t) and particularly dependent on the remaining arguments. 


8. Some Particular Cases 


First we shall consider the case in which we assume in (6.13) that o;,; is inde- 
pendent of Dg,,/Dt, ..., D"g,,/Dt". Accordingly, we take Dg, (Dt=D%g, (JDP = 
--=D"g,,/Di"=0 in (7.7). We see that in this case the constitutive equation 
becomes 
OX; Ox; 
i Aaa ba au mn pa(0 , aah (8.1) 


This result has already been obtained in previous papers [J, 3]. 


Next, we consider the case when, in (6.13), ;; 1s me ordinary function of 
0%,/0X,, Dep f{Dt, ..., Dg, {De, and ae of 6x,/0X,, A Wer wey AM oats arises 
from the eo uaon in (6.1) that o;; is an ordinary function of the gradients 
0% (OX g, US 10%; zeus OD 10 te. We thon have from (7.7) and (7.8) 


oO; bees as i P8pq Aes 

DONO kan UNChe De at. alae (8.2) 
Luendies 200i (¢ ED ORE OL OOS HE ) 
TEES a aN OE EON ANON Sekt Wii Oli, Seix emey ss 


We now consider the further particularization in which F;; in . and (7.2) 
are independent of 0x,/0X,, so that in (8.2) o;; depends on AY ve Ag er ae AG 
only. The particularization can be made formals in (6.1) by taking 0x,/0X;=0;;. 
We therefore introduce this substitution into (8.2) and obtain 


Op lA geen, eel oe (8.3) 


This result has been previously obtained by Riviin & ERICKSEN [6]. Alter- 
natively, we may derive this result by assuming that the stress components o;; 
at time /, in the rectangular Cartesian coordinate system x, depend only on the 
gradients of the velocity, acceleration, ... in the system x. Equation (8.3) may 
then be derived by the method employed in obtaining (5.6) from (5.1). 

The particularization of the constitutive equation to the form (8.3) implies 
that the material is initially isotropic as has, in essence, been pointed out by 
Nott [7]. In [6] and later papers [8,9, 10, 11] further implications of this 
isotropy with regard to the form of the constitutive equation have been explored. 

As a final special case that merits attention, we consider that in (6.12), the 


particular dependence of F;; on 0x,/@X, is of such a nature that the equation 
may be written in the form 
t 
n= Fy ae AD) A aes AS (8.4) 


t=0 
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This constitutive equation can be obtained formally from the more general 
constitutive equation (6.12) by taking 0x,/0X,=6,,. Accordingly, introducing 
this into (7.8) to eliminate 0/0X,, we see that (7.8) becomes 


Crpsewts lotaealles! veel palates 2b (8.5) 
™=0 
where 
Ppa el St) (8.6) 


i Ox» OXy 
We thus conclude that a constitutive equation of the form (8.4) can be written 
in the form (8.5), where F,; is a continuous functional of y,,(r) and an ordinary 
function of AM), ..., Aw. 

This result can, of course, be obtained by assuming that in a specified co- 
ordinate system x, the constitutive equation has the form 


t 


es? OXp(T) | dup” dup” dup” (8.7) 
7 CE MBLOL, A Oh eed, OK, |) ‘ 
t=0 


Then, in the manner adopted in § 5, we find that it must be expressible in the 
form (8.4). We now consider successively deformations which differ from the 
assumed deformation by a rigid rotation at each instant of time except the 
instant ¢. Since the stress components in the system x are unaltered by such 
rotations, we readily obtain the form (8.5) for the constitutive equation. 

The particularization of the constitutive equation to the form (8.5) implies 
that the material is initially isotropic, from considerations similar to those used 
by Nott [7]. The further implications of this isotropy with regard to the form 
of the constitutive equation can be analysed in a manner similar to that adopted 
in [2] in a slightly different context. 


In the case when the dependence on AJ), ..., A\”) is omitted, (8.5) takes 
the form : 
os= FT)  Pya) |: (8.8) 
t=0 


This form of constitutive equation has been particularly discussed by Nott [4]. 
It may also be derived from (8.1) in a manner similar to that used in deriving 


(8.5) from (7.8). 


9. Alternative Formulation of Particular Dependence 


We have seen that if we take as our starting point a constitutive equation 
of the type (6.1), there is no substantial loss of generality in considering that 
the constitutive equation has the form (6.16). Let us assume for simplicity that, 
in (6.16), the dependence of o,, on 0x,/0X,, dv/OX,, ..., dus?/OX, is polynomial 
dependence, the coefficients in the polynomial being continuous functionals of 
dx,(t)/0X,. Then a typical term in the expression for o;; is 


(a) (B) ) 
du, Ou; OU; 


t t 
ig 8S oy, ae Agee ern, ona eRe gg OME s bes e201 BS) ON 
OX;, OX;, OX jn / J : 


Bxp.(t4) Oxp,(t2) Ax (ts) (9.1) 


Cag ero Ng. ax 


UT, OGomonit te) 
qs 
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where the kernels are continuous functions of their arguments, «, f,...,y are 
chosen from the integers 0, 1, 2,..., and v denotes x;, provided the motion 


starts sufficiently smoothly. 
Denoting the Dirac delta function and its derivatives by 6(7), 0’(),..., 6°) (z), 
, we have 
(0) 
a =f t—t) — 2) ee i ae 
OX, 
and (9.2) 
Up 
1 fo" ae (9452.0, 


xX, 


Using these relations it is ETS that we Pe! omit the explicit dependence 
of the functionals on 0x,/0X 4, du NON nec. OX, in (6.16) and on 6x,/0X,, 
avy HOLD Cn. 10%, in (6.17) 1f we allow fae ae to depend on the Dirac 
delta fanction and its derivatives at time ¢. 


10. Relation to a Previous Theory 


RIvLIN & ERICKSEN [6] have investigated the implications of an initial 
assumption that the constitutive equation takes the form of a relation between the 
stress components o;; and the kinematic gradients 6%,/0X,, dv\)/Ox,,..., dv$/Ox,. 
It is easy to show that in a wide variety of cases the assumption of a con- 
stitutive equation of the form (6.17) is equivalent to the assumption made by 
RIVLIN & ERICKSEN. 

We assume in (6.17) that the functional dependence is continuous in the 
range 0S t<¢# and that o;; is particularly dependent on the values of the gradients 
at t=t. We also assume sufficient regularity so that o;; may be uniformly 
approximated by the sum of a number of terms of the form 


ff- x Ghans im iu pla teen) Ti tarts hE 


0 (10.4) 
0 é a é 
¥p, (T%) Xp, (T. ) Xpp (TR) dtd... od the 


Pe he RR Ge 


the coefficients of which are ordinary functions of OH, (OX, COE? [O%q,.. rou HO 
It is easy to show that if the functionals are of the boda type, ihe typical 
term (10.1) must take the form 


t t t 


Sa ATUMET PT Ce ey oe x 
t—T t—T t-T 
x 2%p,(%) Ba 0% pp (TR) at, Ot at irae 
OX, 0X,, OXop 14T,...aTp, 


where T is a time such that y, _ .qn Vanishes if any of the variables 1,,..., tg<t—T. 


Now, if in (2.2) x, is a polynomial function* of 1, it may be expressed as a 
polynomial function of ¢—T, thus: 


%p(t) = %» + (t—z) of) + — (¢—1)20@ +... (10.3) 


* This cannot, of course, be the case if T >. 
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Introducing (10.3) into See 2) it is apparent that (10.2) may be expressed as a 
polynomial in 6x,/6X,, dv\/dx,, ..., dvi*/ax,. Hence, the stress components may 
be expressed as functions of these eeidients 

More generally, if x, may be expressed as a Taylor series*, about the instant 
of time ¢, by an expression of the form (10.3) together with a remainder term, 
and if this series converges throughout the interval t— T<1<t, then we have 
the result that the stress components may be uniformly approximated by functions 
of 0x,/0X,, au) /0x,, a Ov) /0x,. 

If we do not make the assumption that the functionals are of the hereditary 
type, then if x,(r) may be expressed as a Taylor series which is uniformly con- 
vergent in the oe O0StSt, o;; may be uniformly approximated by a function 
of 6x,/8X,, dvP/ax,, ., oven, and t. 

In a similar way we can show that an assumption that the constitutive 
equation takes the form 


t 
B 94,41) (n) 
Oxp,(t) dv ou 
= p ‘ p p 
Bieati| SOP MEL el Us or op os Rally (10.4) 
OX OX OX 


is, in a wide variety of cases, equivalent to an initial assumption that the stress 
components are functions of the kinematic gradients dv!)/x,,..., dv'"/dx,. 
Again, assuming that the dependence of o;; on the kinematic variables is suffi- 
ciently regular and of the hereditary type, we may express o;,; as the sum of a 
number of terms of the type 


t t t 
f fo | tian. -tmyiupih- phan) lat — Ua, +7, tb — TR) x 


EG eC eg ec, avy, (tR) F , (10.5) 
é os ee ae 
OX, OXap 
with coefficients which are functions of CP kaae iB Ove. any. TOs) tas ce: 
time such that 7;..4, vanishes if any of the yasaviee Tae Uo.s ee gail ee T. Again 


we assume, as in the case previously considered, that x,(t) can be expanded as 
a Taylor series about t =#, which converges throughout the interval ¢—T<1z, or 
as a special case can be expressed as a polynomial in T. It is then seen that (10.5) and 
hence o; ; may be uniformly approximated by a function of Gv§)/dx,,..., du}? /6x,. 

CoLEMAN & No. [12] have calculated, on the basis of a emg theory 
the forces necessary to support certain simple states of flow in an isotropic visco- 
elastic fluid. These types of flow had previously been discussed [73] on the 
assumption that the constitutive equation takes the form, advanced by RIVLIN & 
ErIcksEN [6], of a relation between the stress components o;; and the kinematic 
gradients du 0x, ae oul /dx,. It is easy to see that for these flow fields the 
“memory” theory is equivalent to a theory of the type previously advanced by 
RIvLin & ERICKSEN. 

In the case of rectilinear flow through a straight tube of circular cross-section, 
taking the x, axis parallel to the length of the tube, we have 


H(t) =%, %_ (tT) = %, %y(t) = Xs — (¢ — 7) vy”, (10.6) 


x In the case when T>1, this restricts the flow to start at t=O sufficiently 
smoothly. 
Arch. Rational Mech. Anal., Vol. 4 27 
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where v! is the velocity along the tube of a particle of the fluid. v) is, of course, 
independent of t but dependent on x, and x,. We note that in this case pats 
is a polynomial in ¢ —T. 

In the remaining two cases discussed by COLEMAN & Nott [12], Couette flow 
and torsional flow, it is easily seen that, if the x, axis is chosen to coincide with 
the axis of flow, the path of a generic particle of the fluid is described by 

x7 (1) =x, cos @(t— tT) —2, 5m 0? — 7), 

%_(t) = %,sin w(t —t) + x%,cosw(t—T), (10.7) 

By Tie Mears 
w is independent of t and is a function of («{+ 3)? in the case of Couette flow 
and of x3 in the case of torsional flow. We note that the Taylor series expansions 
about t=t for cosw(t—t) and sinw(t—t) are uniformly convergent for all 
finite values of {—T. 

It follows that as far as discussion of these problems is concerned, an initial 
assumption that the constitutive equation for the visco-elastic fluid takes the form 
(10.4) is equivalent to an assumption that o,;; depends on CE eG (Cl 
A similar result obtains in the case of helical flow of the visco-elastic fluid dis- 
cussed in a previous paper [13]. 

Note: In general, equations (6.12), (6.13), (6.16) and (6.17) also depend explicitly 
on the time ¢, except when the functionals in (6.11) and (6.15) are of the invariable 
hereditary type. 

The term hereditary in this and previous papers is restricted to mean hereditary 
which is of the invariable type (the case of the closed cycle). 
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The Effect of an Aligned Magnetic Field 
on Oseen Flow of a Conducting Fluid 


G. S. S. LUDFORD 


Communicated by HILDA GEIRINGER 


Slow flow of a viscous conducting fluid past a sphere is considered. The classical 
Oseen solution is modified by a magnetic field which, at infinity, is uniform and in 
the direction of the free stream. If the pressure of the ambient magnetic field is 
smaller than the dynamical pressure in the free stream, there is no effect on the 
Oseen value of the drag. If it is greater, however, the drag increases with the strength 
of the applied magnetic field and the conductivity of the fluid. 


1. Introduction 


We consider the slow flow of an incompressible, viscous, electrically conducting 
fluid past a sphere, which, for simplicity, is assumed to have the same magnetic 
permeability as the fluid. The applied magnetic field is uniform but weak, 
and is aligned along the free-stream direction. 

This problem has been studied by CHESTER (1957), who neglected the dis- 
turbance of the magnetic field and the inertia of the fluid. This amounts to 
assuming that the Reyneolds number, &, and the magnetic Reynolds number, 
Ry, are small compared to the Hartmann number, /, though the latter is itself 
supposed to be small. In this way he found the drag formula 


D 3 Ii 2 4s 3 4 
(1) “Ds ia 8 i 960 Me 7680 foe eS 


where D;= 6a 09vaU is the Stokes value for the drag in the purely viscous case. 

The purpose of this paper is to take into account, to some extent, the effects 
neglected by CHESTER, by using an Oseen-type approximation in which quadratic 
terms in the disturbance quantities are neglected in the equations of motion. 

As has been found by CARRIER & GREENSPAN (1959) for a plane problem, 
the ratio “ H/o) U2= M?/RRy of magnetic pressure to dynamic pressure is an 
important parameter in the discussion. For M?<RR,,, the two wakes which 
appear both lie downstream of the sphere, while, for M*>RA&),, one les up- 
stream; and the drag experienced by the sphere depends in a completely different 
way on M, R, and Ry, in the two cases. Indeed we find 


D 3 
(2) De — 1 + i I, 


Di 
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where ie for M*<RERy, 
a 2MPt BEM fog MARR 


and K is correct, in the present approximation, to the first order in M, R, Ry*. 
Higher-order terms are irrelevant, since it is known [PRoUDMAN & PEARSON 
(1957)] that, in the purely viscous case (M=R,=0), the Oseen theory only 
gives D/Dg correct to the first order. 

For M=Ry=0, K=R and we obtain OsEEN’s result. When R=Ry=0, 
K=M< and (2) gives the first two terms in CHESTER’S formula (1); the remaining 
terms correspond to higher-order terms in K, and are therefore suspect for the 
reason just given. 

The surprising thing about (2) is that the magnetic field has no first-order 
effect on the drag so long as M?<RRy, but that for M?>RR,, it does. In the 
latter case K is always greater than R: as might be expected, the attempt of the 
magnetic field to inhibit motion across its lines of force, as the fluid passes around 
the sphere, results in an increase of drag. 


Similar results have been found by VAN BLERKOM (1959). 


2. The Equations of Motion 


The steady motion of an incompressible, viscous, electrically conducting fluid 
of constant properties is governed by the equations 


3) (a) v-gradv=—— gradp+V?0+ © curl xH, (b) divw=0, 
0 0 
(c) curlH=o(E+ypux), (d) cul—=0, (e) dvH=0. 


For an axially symmetric motion in which the velocity v and magnetic intensity H 
le in the meridian plane and are independent of the azimuthal angle, the con- 
duction equation (3c) shows that the electric field E is perpendicular to this 
plane and also independent of the azimuthal angle. Equation (3d) then shows 
that H=0, if it is to be finite at the axis. 

When v, r, and H are now made dimensionless by referring them to the 
velocity at infinity U, the radius of the sphere a, and the external field intensity 
H,, the equations reduce to 
(4) (a) Rv-grady =— gradp+ V20+ Re colHxH, (b). divv=0; 
(c) culH=RK,voxH, (d) divH=0. 


Here R= Ualy is the Reynolds number, Ry=Uapo is the magnetic Reynolds 
number, and M=wH,a(o/g9v)* is the Hartmann number; the pressure is now 
given by (@9v U/a) p. 

If the x-axis is taken along the axis of symmetry, then at infinity v= H=i. 
Following OSEEN’s work in the purely viscous case, we linearize these equations 
by writing v+t for v and h+ié for H, and neglecting quadratic terms in v 


* In fact, this result is also valid when the permeability of the sphere is different 
from that of the fluid. 
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and h. From (4a) we then find 


yp ain a ik REE, 
(5a) Viv —R 4 oe ae (p pate i), 
while the curl of (4c) yields 
2 ow Oh\ _ 
(5b) Vh + Ry (2° — "| =0, 


when (4b) and (4d) are used. This pair of equations may be replaced by an 
equivalent pair of Oseen equations, so that their solution follows a well-established 
path. 


3. Solution of the Oseen Equations 
From (5a) we have 


2 
V3 (b+ zhi) =; 
and a particular integral of equations (5) is therefore obtained if we write 
Ry? + M?h-i=(M?— RRy) sep v=h=gradq, 
where @ satisfies 
(6) V2p=0. 
The solution is completed if we set 
w=gradg+v’, h=gradgp+h', 


where wv’, h’ are solutions of the corresponding homogeneous equations. Two 
linear combinations of the latter are the Oseen equations 


O- 
2 Ox 


- o- 
1 Ox 


(7? 2k )(e' on == 0% (v? 2k J (e’ | ¢ hh’) =0, 


where 


hy, ha == [R+ Ry F(R — Ry)? + 4), 


Hy» Ay = ae [AR Rye /(R Ry)? 4 4M? ); 


in addition we clearly have 


div (v’ + a,h’) = div(w’ + a,h') =0. 
It follows that we may write [LAMB (1911, 1932) | 


v' + oh’ = (%) — 04) 


1 3 
2h, grady,— % i ; 


v' + ah’ = (a, — a) 


1% fille: 
Fee rai, 


(8) (02 — 2h 2) m= (02 — 2h Ge) ta = 9. 
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This leads to the complete solution 


v = grad (p+ a 7 ae ae 1) — (2 %1 — % Ho) 4, 


(9) h= grad (p— oh Near an is| te AYs 5) 


em ROY yaaa ek Sia.) , | 
7" a Ry (es Ox 2h, Ox (a — Xa)]. 


Inside the sphere curlh=divh=0, and hence we may set 


(10) h = grady, 
where 
(11) V2y~= 0: 


The fundamental solutions of (8), from which the general solutions may be 
constructed (Section 4), are exp[—|hk,|r+,x]/r and exp[—|k,|7+,x]/r, 
respectively. These lead to wakes in the directions of increasing k,x and k,x, 
respectively. Now kh, is essentially positive, so that there is always a downstream 
wake. On the other hand, we have k,20 according as M*7SRR,,. Thus, when 
Hy is less than @, U2, the second wake is also downstream, but when it is greater 
the wake is upstream. 


4. Boundary Conditions at the Sphere 


Since, for simplicity, we assume that the permeability of the sphere is the 
same as that of the fluid, the magnetic intensity h must be continuous at the 
surface of the sphere. On the other hand, the disturbance velocity wv of the fluid 
assumes the value —i there. From equations (9) and (10) we see that this places 
the conditions 


Op tan@ eis Aenea _ oy 
or 2h, Olesen (a — Xa) Sos 
oy 1 OX, iL tale ; _ oy 
ah Bh 06 ) Bho Bie oe! eee 
and 
Op OX) % Ox» - = 
or QR ey 2k, Or (241 — % Ya) cosd = — cos®, 


OP Sor oN & Ox, Ax : mye 
a6 ah, a8 a ot eX % %2) sind =sin#, 


on the functions @, 7, 72, y at the sphere. 


Formulas given by GOLDSTEIN (1929) for the purely viscous case enable us 
to obtain the solution of the present problem to any desired accuracy. However, 
we are only interested in the first-order effects in M, R, and R,,. For this purpose 
it is better to adopt the ad hoc method used by Lams (1911) for the purely viscous 
case and later by CHESTER (1957) for the present problem. Its justification lies 
in the fact that, as far as it goes, it is equivalent to GoLDSTEIN’s method. 


Now the general solutions of (6) and (8) outside the sphere are, in the case 
of axial symmetry, 


es) A or 4 woe eo B gn ane b ee) 6” fener 

=> ee = 1 ES = a, —————— 

u Zs ” Oxn ee ey 2 oy” ( Tp ee 2, Ow aan ( 
= n= 


if 
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where the + sign in y, corresponds to k, 20 (the solution must vanish at infinity). 
The general solution of (11) inside the sphere is 


(oo) nee Somer 
= n 
p— > Dit ex" (+) 
n=1 


for the case of axial symmetry. When these are substituted into the four boundary 
conditions and all functions are expanded in powers of k, and ky, and then re- 
arranged into spherical harmonics, we find that A,, B,, C,,, and D,, must satisfy 
the equations 


— Ag Sp — 24 S++ Bo op So+(—1+4 hy) Sy + 4 hy Sp +0 (8) + --- 


— Col sp Sot (= 14 Sha) Si + 4 hy Sp +0(88)| — + =D, Sb, 


Ay Sp+-+- + Bo|(— Stith) S++ hy Sy + 0(88)| +o 


—€,|(— a taza) Sots hy S, + 0(8)] — Se was 


— Ay Sy — 24; Sy ++ = 03 Bo 3 So+(—15-f) Si +m Sp +0(88)|—--- 


+a Co| : 


apy Sot (145 Fa) Sit he So +0 (88) +2 = Sy, 


Ay Spt + ~ ay Bo|(— 5 £ i) Sot > hy Si +O (8i)| — ~~ 


1 
4 


+a Co|(— 5 + 5 ha] Sot 5 hx Si +.0(88)| += Sp. 


Here Sy, S,, S. are the spherical harmonics 


ba yal CZ 
Ox i Ox? Bie : 
and the dots stand for omitted terms in A,, D,, (1>1) and B,, C, (n>0). 


Since the spherical harmonics are independent, these provide an infinite set 
of equations for the coefficients A,, B,, C,, and D,. From this set the A, 
and D,, can be eliminated, the first three of the resulting equations being 


’ 
r=1 


(1+ as) Bae + 0(ii)| a tea) Col ay + 0(%8)| BR RA 
og By| 2+ 3h + O(AH)| +++» — 01 Co] — 2+ 3 he +0(88)] — = —3, 
(1+ a) By| 4 + O(A8)] + ---— (+04) Co] > +.0(88)| — = 0. 


Sufficient accuracy is obtained by solving the first two of these equations with 
io O = OL0t hn 0: 


B=Q=>— 


Jae Oy ky — a Re) | i 
4 (%— 0%) 
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The bracket is correct to the first order in k, and ky. The corresponding value 
of A,, which, as we shall show in the next section, determines the drag, is 


yee 3 (ky — ho) | 3( 4%, moet | 
12 2 4 (2 — Oy) Ry Ry 4 (a2 —%) 
ve) ks oi [1+ 2K 
2(R Ry— M7?) 8 , 
where Reviiors kppaoy 


K=4 2M*+R*—RRy 
/(R—Ry)?+4M? 


for k, <0. 


5. The Drag on the Sphere 


The drag D due to the pressure and viscous forces on the sphere is @) U?.a?/R 
times 
(13) —fllptmi—nnjds, 
[cf. LamB (1932)] where (J, m, 2) are the direction-cosines of the outward-drawn 
normal to the sphere S: y=1, while 7 and € are the y- and z-components of 
curlw. The drag due to the Maxwell stresses is zero: it is proportional to the 
flux of h through the surface of the sphere. 


Now, according to equations (9), the integrand in (13) can be written in 
the form 3 Me 
Op 1 Oe) f 90g 
Altes shag oa ers De) ak a va w2)} | 


Ox OX Ox OX 
| mm i 2 t. he 2 
jm | ap ae | te [oa 


while on the sphere equations (7) show that 


LD) 


1 : 1 , - 
per l= eee aad ah-+ (a+ Ry) grad ¢| 


for k=k,, k, etc. Hence the integrand is equal to 
_ M*\ ap M2 
- [R Raa) Ove Ry hy 


where /, is the component of h along the outward normal. Clearly the h,-term 
makes no contribution to the drag, and 


ed We MALE R Me Spise avon M2 
Ds 6x0,U2a/R rele Rhee earl RRy) 40 


where D;=6x@)vaU is the Stokes drag. Inserting the approximation (42) 
for Ay, this becomes D 


He = 1+2%R, 
ute | 1 for M®<RRy, 
ee) ree ————— leiorh Mee 
| R LV(RiRy —VRy (RY + 4 MR Ry i 


Thus, when the magnetic pressure ju H@ is less than the dynamic pressure 0) U? 
(1.e. M?/RRy<1), the magnetic field has no influence on the drag. The effect 
when “ H¢> 6, U? is illustrated in the figure, where x has been plotted against 
M*/R Ry, for various values of Ry/R=o. 
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In the purely viscous case, o=0, M=Ry=0, and x=1, so that we recapture 
OSEEN’s result. On the other hand, when R and Ry are negligible, x= M/R 
and we obtain the first two terms of CHESTER’s (1957) formula (1). In the figure 
this corresponds to the asymptotic behavior of the curves for M2/RR,, large. 


4 = 
15 |- LOV=09 
L YOUY9 
7a 
L 
70\— 16 
|= 
| 40/43 
L y 
a 16/9 
7 
| ey 
q é LEE TA 
7 - =0 
eee eee ee | ! ! | | 


5 0 
wH3/QoU" 


Fig. 1. Variation of the drag factor x with HHS /0o U? for various values of uo » 


Note added in proof. As was stated in the introduction to the paper, Dr. R. van 
BLEeRKom obtained similar results in his Ph. D. thesis at Harvard University. This 
was written under Professor G. F. CARRIER and is dated May 1959, some two months 
prior to my starting work in this area. In the same month, Professor CARRIER very 
briefly presented the results at a meeting I attended. 

I had been made aware of these facts just before submitting the paper to this 
journal. At that time it appeared that Dr. VAN BLERKOm’s interests were being 
safeguarded by his concurrent publication in another journal; there was no doubt in 
my mind that Professor CaRRIER’S announcement had made no impression on me. 
However, since VAN BLERKOm’s publication has been delayed, it now seems appro- 
priate to make the situation clear. 

Research sponsored by the Office of Ordnance Research, U.S.Army, under Contract 
DA-36-034-ORD-1486. 
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The Method 
of Extremal Points and Dirichlet’s Problem 
in the Space of Two Complex Variables 


JERZY GORSKI 


Communicated by S. BERGMAN 


Introduction 


NEVANLINNA’S theory yields information about the distribution of values of 
a meromorphic function of one complex variable. In an attempt to generalize 
this theory to functions of two or more complex variables we meet two essential 
difficulties: 1° the zeros of analytic functions of two complex variables he on 
segments of analytic surfaces, not at isolated points, and 2° the formulation 
and the solution of a boundary value problem analogous to the Dirichlet problem 
is not immediate. In contrast to the case of one variable there does not exist 
in general a pluriharmonic function (¢.e. the real part of an analytic function of 
two complex variables) which assumes prescribed values on the three-dimensional 
boundary of the given domain. In connection with the second problem S. BERG- 
MAN classifies the domains in the space of two complex variables (see [/]) into 
4° domains such that to every point P of the boundary there exists an analytic 
function which assumes at P a value greater than that at any other point (e.g. 
the hypersphere), and 2° domains where only a certain part of the boundary 
has this property; it is called the maximum boundary manifold and is denoted 
in the literature as the Bergman-Silov boundary. Further, BERGMAN showed 
that there exists a class of domains for which the maximum boundary manifold 
F is two-dimensional?! and to every real continuous function given on F and 
satisfying certain additional conditions there exists a pluriharmonic function 
which is uniquely determined in the four-dimensional domain considered and 
which assumes on J’ the prescribed values (see [2]). 

Since not every given boundary function satisfies the above mentioned 
additional conditions, BERGMAN introduces a larger class of functions (functions 
of extended class, see [3]) than the pluriharmonic functions; use of this class 
enables solution of the boundary problem with values prescribed on the Bergman- 
Silov boundary. Introduction of the functions of extended class can be effected 
in various ways. In connection with these investigations D. B. LowpDENSLAGER 
[19] and H. J. BREMERMANN [4] defined different extended classes of functions. 


‘In this case F is called the distinguished boundary surface. 
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The aim of the present paper is to show how we can obtain a certain extended 
class of functions using the method of extremal points introduced and developed 
by F. Lrya (see [13], [14], [15]) and others ([7], [8], [9]). The class obtained 
is a certain subclass of plurisubharmonic functions, and the solution of the 
boundary-value problem considered is expressed as a generalized potential of 
a simple layer (see p. 418). 

In Section 1 we generalize the method of extremal points to the space of 
two complex variables. In Section 2 a specialization of functions of extended 
class is made. In this latter case the original proofs based on the use of Stieltjes 
integrals can be replaced by an elementary method. An example illustrating this 
approach is given. 

It should be mentioned that our results can be generalized to the case of 
analytic functions of m complex variables. 


1. The method of extremal points in the space of two complex variables 


As indicated in the introduction we use the method of extremal points in 
order to construct the extended class of functions. First we introduce certain 
notions; later they will be specialized to the case in which we are interested. 


1. Let E be a compact set in m-dimensional Euclidean space R,, and let g(, q) 
be a real, continuous function of two points # and g defined in R,,. Suppose 
that g(p, 7)=¢(¢, p) 20 and g(p, f)=0. Let f(f) be a real, continuous function 
defined on E, A>0O a real parameter and pM={p™ po”, ..., o£) an arbitrary 
system of »+1 points of E. A system q™ (A) = {q, gq”, ..., g?} of n+1 points 
of E is called an extremal system if the inequality 

V,(9(A)) = IT 8 (a;”, a) exp {— Alta” ee) 
eee 
holds for every system p™€ E. 
It is known (see [/3]) that the limit 
2 
lim [V,. (a (A)) | WOT i= (i =o 0A) 


n> 


INV 


0) 


exists, and we shall call it the generalized capacity of the set E. 

Special cases. 1. When E is a plane set and g(p, g)=|z—¢], where z and ¢ 
are two points in the complex plane, then d(E,|z—€|, 0) is the transfinite 
diameter of the set E (see [5]). 


2.In the case when E=R, and g(f, q) =exp {— Al where fq denotes the 


Euclidean distance between two points # and 4, aE, exp ewan 0) is equal 
to exp {- <i} where C (E) is the Newtonian capacity of the set E (see [2/]). 

3. In the case when E is the space of two complex variables z, w and g(p, ¢) = 
1\zy —wl|,d(E,4|zy — wl |, 0) iscalled the triangle capacity (1’écart de l'ensemble) 
of E (see [16}). 

The first two kinds of capacity play an important role in the theory of con- 
formal mapping and D1rIcHLET’s problem; the last one has been used to de- 
termine the domain of convergence of a series of homogeneous polynomials. 
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The inequality d(E, g,4)>0 depends on EF and g(f,q) but not on 1D 
d(E, g,0)>0, then d(E, g, A) is positive. Suppose that the set E and the function 
g are such that d(E, g, 0)>0. 

To each system of extremal points g(A) there corresponds a function of 
the Borel set e defined by the relation 


0, if e does not contain any point of g (A), 
(g) — a 
Haale) = Maal?) _! if e contains k points of g™ (A). 
n+ 
The functions y,,,(e) satisfy the inequalities O<w,,(e) <1 and therefore define 
a certain distribution of unit mass on E. Let mw, be the limit of a convergent 
subsequence chosen from {w,,,}; then 


(1) logd(E, g, A) a, log g(p, 9) du, (P) du; (9) — 24 J f(a) du, (q) =I(u,) 21 (z), 


where 7 denotes an arbitrary distribution of unit mass on E [see [7], where the 


special case g(p, g)=exp {- mr is considered) : 
Consider the function 


v,() mail log g(P, 9) duu, (q) — Af (A) 


defined on E, and let E,(g) =F, be the support of the mass distribution yu‘! = m,. 


Theorem 1. The function v,(p) 1s equal to a constant almost everywhere on E, 
(see [6]). 

Proof. Set? 

A v,(p) dua (b) = 42. 

We cannot have v,(p)=y,-+¢, ¢>0 at every point of EF, and hence there exists 
a point p,€ E, such that v,(9)<y,+¢e. The function v,(p) is upper semicon- 
tinuous, and therefore this inequality holds in a neighborhood O (A,) of the point fy. 
Suppose there exists a set F CE such that d(F, g, 0)>0 and v,(p)=y,+2e on F. 
Set 4, [0 (fo)|=m, and let o=a(e) be a function defined by the formulas 


a(e)=— (2) on O(p), 
a(e)=0 on Ff, 
o(F) =m 
a(e) =0 outside O(f)) + F 
and 
A= 


a log g(p, g) do(p) da(q) > — &. 


E 


* Almost everywhere means except for a set Fc E, such that CE TINO)! 10), 
* va = log d(E, g, 4) + 4 J £(4) dy (q). 
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The function »=y,+ho, 0<h<1 defines a new distribution of unit mass on E. 
From (1) we obtain 


feat) log g(p, g) dv (p) dv(q) — 24 J f(g) ae(q) 5 
— J floge q) 4, (Pp) ding )-+28S 9 Jaya, ( (7) 
=hA + any v,(¢) da(q) 


Z2hA +2h[—m(y, +6) +m(y, + 26)] 
=h?A+2mhe. 


For sufficiently small h>0 the right side is > 0, in contradiction to the inequality 
on the left side. 
Therefore v,(p) <y, almost everywhere on E (cf. footnote 1). We cannot have 
0, (Po) <<, at a point f)€ E,, since because of the upper semicontinuity of v,(f) 
this inequality would hold in a neighborhood of #, and f v,(p)du,(p) would be 
E 


<y,. Hence v,(p)=y, almost everywhere on F,. Q.E.D. 
Let u,(p) be defined as follows: 


0(6) = | loge, 9) d(a). 


From the upper semicontinuity of u,(f) it follows that 
(2) lim U,(p) Su, (po) = Y2 1 #(p)) almost everywhere on E,) ‘ 
P— po CE, A 


Suppose g(p, g) +0 for +g and log[1/g(f, g)] satisfies the maximum prin- 
ciple of O. FRosTMAN, 1.¢. that for every ese o of positive mass on £, from 


the inequality of 
+ flog do(g)SK<~% 


on the support £, results the same inequality in the whole space. Then 


(3) lim = ftoge(.9) ding) =? + Ibo). 
E 


p> br C Ey 


For let #, be a point of E, and s, S be two spheres centered at p,) such that 
0> 4 | loge(po,9) du,(q)>—é 
s 


and 


7} eel ’, 9) diy (q ~ 7 {sel P'.9) diag a) <e 


for p, p’Es. Theretore we have the wee inequalities: 
4 1 Sac el 
(i) + [ osete, q) dunia ,ink Ly frees. 9) an} > 


> flosett.9) du,(q) —€ 
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where #, is a suitably chosen point on £; s; 


(ii) ne al log g(P, 9) dala al log g (Po, 9) 4ua(g) > 


> § fss (Po, 9) oalg) > —€, 


(ii) flog g(p,9) dua(a) > x} tose pr,9)dyn(g)—8, BES. 


B=Ss 3S, 
(i), (ii) and (iii) yield 


lim 4 f log e(,9) dia(9) zs | ssl pr,9) dy14(q) — 36. 
E 


p— po CE, 
Therefore 
lim + flog g (2.4) dya(a) = lim > flog g(61,9) dala) = 7 +1 (bo) 
p>poCE, A — Po d. 
E pee; E 


almost everywhere on F,. 
From (2) and (3) we deduce the existence of the limit 


lim , 4 flog ¢(b, @) dra(a) = 7% + F(bo) 
E 


p> Po CE, 


almost everywhere on £,. 


When the function log [1/g(f, q)] satisfies the principle of energy (see [6], 
[20], [22]), that is if 


| fre AG 
Eads 


is always =0 and is equal to 0 only when t =0, a 1° the function uw, is unique 
(and therefore the sequence {/,,,} is convergent), 2° 4,49 when A—>0 (see [12]) 
and therefore F,—>E,. 


1°. In fact, if two different functions o, and mw, realised the upper bound of 
the oe 


t(p)dt(q),  t(E)=0, 


= J floss (p,q) dt(p) d q— 24 J ta) 
then 


esti! log ¢(p, g) do, (p) do, (q )— 24 fila) ) do,(q) — 
~ f Jog ¢(P, 4) dua(b) dynla) + 22 J Hla) din(a): 
Set o,=/4,+». We obtain 
O=— 24 J t(g) dv(q) veal! log ¢ (p,q) dv(p) dv (gq) + a) log g(, g) du, (Pp) 4v(q) 
= a vag) dv(q) + J log g (b. 9) dv (p) dv(q) 
S2[y,»(E,) +y,9(E — E,)] lal log g(P, q) dv (p) dv(q), 


Functions of Two Complex Variables 417 


and, since »(£,)-+»(£ — E,)=0, it follows that 
Dea) log g(p, q) dv (p) dv(q). 


log |1/g(P, ¢)] satisfies the principle of energy, and therefore the right-hand side 
of this inequality is <0, the sign = holding only in the case y=0. Therefore 
Oo, =, - 

2°. Suppose /,, 0. We have 


logd(E, g, 4,) —logd(E, g, 0) =v a log g (p, g) tt, (b) 4p, (7) — 


— dn J Ha) dita, (g) — J Jog (. 9) 4H0(b) di0(9) 
Set ¢,=/M,, — My; then ‘ 


logd(E, g, 4,) —logd(E, g,0) =f flogg(p,¢)do,do, + 
EE 


ae 2 i. ; log g (2, q) dig do, — 2d, J £9) Aly, 


The left-hand side tends towards 0, when -> co, because d(E, g, A,,) is continuous 
with respect to A, (see [11], [12]). On the other hand 


aj J logs (2, q) dig do, S 2logd(E, §) 0): [o,, (Eo) dsOn\ —L5)| = 0 


Therefore 
0< lim fat log g (p,q) do,,do,,. 


N—> OR 


Let o be the limit of a convergent subsequence chosen from o,; then 
OS J floge(p,9)dodc. 
EE 


From the principle of energy it follows that equality can hold only when o=0. 
Remarks. A. For some special cases, e.g. n=2, g(p,q)=|z—¢| or n=3, 
g(p, 9) =exp {— ahs the continuity of the function u,(p) has been proved when 
E is sufficiently regular (see [/4], [10]). Similarly in these two cases it has been 


proved that under some additional conditions on the set / and the function /() 
there exists a A)>0 such that E,=E for A>0, A<A, (see [9], [10)). 


B. Let E be the boundary of a bounded domain DCR,,, and let 


p(w) = D aselP) ult 0 (6) +S ba(B) w, (6) +0(8) u(f) =0 


be a linear differential equation of elliptic type. Suppose I’(#, ¢) is the fundamen- 
tal dee of y(w) =0 which satisfies the following condition: I"(f, 7) =I"(¢, f) > 0, 
= On) Set 217.0) )=exp{— I'(p, ¢ )} and suppose I'(p, g) satisfies the 


re Pp) 
maximum principle of peered: If d(E, exp{— I'(p, q)}, 0) + 0, then the corre- 


sponding function u,(p) = f I’(p, g) du, (gq) is the solution of the equation p(u) =0 
E 


which assumes on £, the value . + f(P). 
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C. Let D be a four-dimensional bounded domain in the space of two complex 
variables z, w, and let /(p) be a real continuous function defined on the boundary 


Biot DT Set 
g(b, 9) =|h(b,9)| 

where h(p, q) is a continuous function of two points p=(z, w), g=(¢, 7) defined 
in a domain! G>D+E and satisfying the following conditions: 

1. |h(p,9)| =|h(g, )| 290, 4(b, 2) =0, 

2. h(p, q) is an analytic function of p=(z, w) for fixed g=(C, 7). 

Suppose that d(E, ||, 0) is positive. Using the method of extremal points 
indicated above, we obtain the function 


; 
u(6) = [ log] A(,9)| des(a) 
E 
which is plurisubharmonic? at every finite point. If h(p, g) +0 for Prey and 
q€E£, then u,(f) is a pluriharmonic function in D. We have 


| < 4 (p) almost everywhere on E, 


= 424 #(p) almost everywhere on E,. 


If the function log c : satisfies the principle of energy, then 4,» for A—0 


1(p,9)| 
and therefore E,->E, (see [12)). 


From the upper semicontinuity of the function u,(f) follows 


©) : tim pj talP) S a (Po) 


We shall prove that for every “‘regular”’ point py € E, the limit 
lim (2) = (Po) 


p> bo CE, 
exists. 
Let py be a point of E, and let F(fp) = E> {h(f, 7) =0} be the product of 
the sets E and {h( (fo. 9) =0}, From the aes 4 the function h(f, q) it 


follows that for | pp 9|<7, the set En{h(p, g) =0} is contained in the set G;(po) = 
EOE, (fo), where Fy (py) denotes the hull of the set F(f)) with radius 6>0. 

Let e>0 be an arbitrary number, and let ~, be a point of E,. If dis sufficiently 
small, then 


1 
ay [ 10g] (bo, 4)| diea(a) <6. 
Therefore ss 
(ES ine il log |h(p, q)|dua(q = log |h (Po, g)| du,(¢) > Ua (Po) — 
pb. CE; - al Z § 0 FZ) FAW a(Po 
E—G; (by) E—G5 (Po) 


'G can be the whole space. 
a A real function V(p) is plurisubharmonic in a domain if the following conditions 
are satisfied: (a) — co <V(p) < ce, (b) V(p) is upper semicontinuous, (c) the restriction 


of V(p) to any analytic plane P of complex dimension one is subharmonic in the 
intersection PD. 
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We shall call a point )€ E, “‘regular”’ if for every ¢ > 0 there exist two numbers 
6, (Po) > 0 and 7(e) > 0 such that 


(#8) xf elt aldrala|<e for |p Po] <r(6). 
G 6¢ (Po) 
If PoC E, is regular, then from (**) and (***) follows 
lim (Pp) 2, (bo) 5 
P— bo CE, 
and therefore (see (*)) the limit 


lim _ 1,(b) =u, (Po) 


P— Po CE, 


exists’ -OrE7D: 


Suppose fy€ F, is a regular point with respect to the Dirichlet problem for 
the domain D, and suppose there exists a cone c contained in D with the vertex 


at py. We shall show that u,(po) = a + (po). Let s be a sphere with radius 7 
centered at p). We denote by u3(p) and wu; (p) the following functions: 


wl) => f loglr(b.g)ldualg), Wi) =+ f 108 |h(P, a) | dul). 
G5 (Po) E—G5(ho) 


Let k be equal to the ratio colume oles 


. Suppose the number 6>0 is such that 


volume of s 
(a) O>u,(by)) >—ke, e>0. 

We choose the radius 7 so small that 

(b) uj(p) >u,(bo) —€ and ,(p)<O for pés. 


Let w(p) be the generalized solution of the Dirichlet problem for the domain 
D with the boundary value /(p). The function ,(f)—w(f) is subharmonic in 


/ 


D, <y,/4 almost everywhere on E. Let m,,, m;, and m,, be the average value of 
u,(p) —w(p), uz (p) and wj(p) on cs respectively. We have (w(p)<f(o) + ¢ in cs) 


4 a Mes = Mes Ms I (Po) é. 
A 


From (a) and (b) follows 


Mi, = Uy (Po) — €- 
Therefore 
PA > ui; Pe) — H(bo) —3 8 2 abo) — I (bo) — 3°, 


and since ¢>0 is arbitrarily small, we obtain u, (9) S 5 +} (po). 


On the other hand we have ty (bo) 2 +1 (bo) for ,€E,, and therefore 


u, (Po) = TA + t (Po) - O° ED: 
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From the definition of the extremal points and from the analyticity of the 
function h(p, g) with respect to =(z, w), it follows that in the case A=0 we can 
choose the extremal points on the Bergman- Silov boundary B of the domain D. 
(In the case of a domain with the distinguished boundary surface / we have 
F=B.) We denote by Ey the support of the mass distribution 7) defined by 


the extremal points in the case A=0. In general Ey +B. If log ae satisfies 
the principle of energy, we have E,—>E, for A—0. [h(, 9) 


We set A=1, and we consider the functions 


ul (p) — 1 = (8) = 7 


which correspond to different functions /(f,qg) satisfying Conditions 1. and 
2. (p. 418). For #€D we define the following function: 


,(p) = Tim {sup [uf (@) —"]}. 
IP \ (nh) 
Since every function u{")(q) is plurisubharmonic in D, the function w() is pluri- 
subharmonic in D (see [4]). 
Suppose the domain D is regular at a point ~),€E (see p. 419). For 
i ae h(p, g) we have wu) (gq) —y<f(q) for EE. If! gc E”, then 
(gy) —y? =f (go) and therefore sup [w\ Nig ‘) —y\"]}=f (qo). This formula is valid 
a every point g,€ >) E™. (%) 


(h) 

Let w(f) be the solution of the Dirichlet problem for the domain D and the 
boundary values {(f). We have u)(¢) —y{ —w(q)S0 on E, and, since the function 
u (q) —y? —w(q) is cihhalenenis in D, this inequality holds in D. Therefore 
ae [{” (gq) —y{? —@ (q)]S0 for g€D and w,(p) —w(f) SO in D, and consequently 


) 
w(p)Sa@(p), PED. 
Let q be a point of the set ae then 


lim w,(p) < lim w(p) =f (%) 


A f ee « PQ P—>% 
and, if gp is “‘regular’’, 


Jen — tn (ap len a tte 


Therefore? for regular go 


@ (P) = F(4o)- 


lim 
Pq cz gE) 


1 EW is the support of the mass distribution gle 


? Analogously we can set A= — 1 and introduce the function 
w-4(P) = lim fint (ol!) (q) — yh} 
q—>P 


which is plurisuperharmonic in D. We have 


w(p)2o(p) and jim w, (b) =f(%), for a “regular” point Woe DED 
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Remark. Let D=D,xD, be a product domain, where D, and D, are simply 
connected domains lying on the z and w plane respectively. Suppose the boundary 
FE, and E, of the domains D, and D, are regular Jordan curves such that the 
angle between the tangent and the real axis expressed as the function of the 
length satisfies the Lipschitz condition. Let /(p) be equal to /,(z)+/,(w), where 
f(z) and f,(w) satisfy the Lipschitz condition, and let h(p, q) =(z—) (w—7). 
Then for sufficiently small A>0 the function w,(f) is pluriharmonic in D, con- 
tinuous on D+F, E=E, XE, and equal to /(p) at every point pCE (see [10]). 


2. Specialization of the function h(p, q) 


The considerations in this section give an explicit construction of a function 
of the extended class (a plurisubharmonic function). 

Let D be a four dimensional bounded domain in the space C, of two complex 
variables z, w and f(z, w) a real, continuous function defined on the boundary 
E of D. Set h(p, q)=2zn—wl, p=(z, w), g=(6,m). It is known (see [16]) that 
d(E,|zy—wl|,0) is >0 if E does not contain the origin (0, 0) of the coordinate 
system and contains two points f=(z, w), g=(¢,7) whose “triangle distance” 
|pq|=|2zy—wl| is >0. We assume that the origin (0,0) lies outside of E. 


We denote by ®” (p; £™), 7=0, 1, ..., m the following functions: 


(4) pp) = 7 Wks a 2 exp [Af (p,)] 
oe 


where pp,=(zw,—2,@). (pb; p™) are homogeneous polynomials of the 
variables z and w. 


For a fixed point p € Cy we set 
(5) Dip: f) ==) int {max |}? (b Os Wi 1, Qe 
(7) 


POCE 
It can be proved (see [16], where the case A=0 is considered) that the functions 
(5) satisfy the inequalities 
D(bi A) 2P,(b; A) P, (254), w= 1,2,.-. 
and, therefore (since d(E, |zy7 —w¢|, 0) >0), the limit 


(6) P(p; A) = lim V®,,(; A) 


exists at every finite point p= (z, w) of C,. 

Consider the functions DY) (p; g”)(A)) where g)(A) denotes the 2" extremal 
system of points with respect to the function h(p, g)=z7—w¢. From the defini- 
tion of the extremal points g(A) it follows that 


(7) |B! (p; g@ (A))| Sexp[nAf(f)],  7=0,4,...,0 
on £. We set 
A (ge(A)) = IT (2; Wp — 2, @;,) exp {— Alf (z;, @) + f(z» a) I}, 7=0,1,...,” 
kj 
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and suppose that 
|A (g(A))] SA (g@(A))| for 7=0,1,-..,” 


At every point (z, w) = of C, which satisfies the condition 


(8) 0(p, E) AS Nil —wl|>0, 
the limit 
(9) jim n |p (pq (A))| = Op; A) 


exists. On the other hand at every point p€ C, the sequence 


{[ 2 eee: anl]t} = foul: a} 
is convergent to ®(p; A) (see [16]): 
(10) lim , (p; A) = O(p; A). 


n— CO 


The function ®(p; A) has the following properties ?: 
4. For every point p=+-0, ®(p; A)>0 and D/P; A) satisfy the inequalities 


Lele lel om < 61; 2) Sse aR 


where m = inf /(p), M=sup/(p) and A=sup (|z], |w]). 
PEE PEE PEE 
2. B(p; A) =0 at f=0 and is continuous there. 
3. For an arbitrary complex number c 
P(cp; A) =|c| D(p; A) 


4. Let S be a set of points of C, which have positive triangle distance (8) 
from the set E. S is an open set because 6(f, £) is continuous function of the 
point ~. The convergence (9) is uniform in every closed set GCS, and therefore 
@(p; A) is equal in S to the absolute value of an analytic function. 


5. The function ®(f; A) is lower semicontinuous everywhere. 
6. On the set E 


(14) O(p; A) Sexp[Af(A)]. 


7. When exp /(p) = 
nomial of degree k, then 


®(p; 7) =exp|Fi(p)| on E. 
8. Let £ be contained in the set of points defined by |Q,(f)|—=c>0. Then 
D(p; A) Sexp[Am] on E. 
9. The function ®(p; A) satisfies everywhere the inequalities 


O(p; 0) exp[Am] S O(p; A) < O(p; 0) exp[AM]. 


A (p)| where Q,(f) is an arbitrary homogeneous poly- 


' Properties 1—9 can be easily proved; cf. [16]. The main tool is the Lagrange 
interpolation formula. 
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Remark. Using the earlier representation of Section 1, we can write 
log ®(p; 0) = J log jan — we | duug(C, 7) — log a(E, |2q — w6 |, 0), 
log P(p; A) ai) log |z — we| dy (6,0) —y, = (h) — Y2- 
Let Ej be the set of points of accumulation of the set of all extremal point systems 
qo, 9, 
qo, gy, 9 


and suppose that every point py)=(%, Wp) CE possesses the following property: 
(H) for every sphere K (po, 7), the set EK (po, 7) contains an arcy of a curve 
CCE, pfyo€y whose triangle capacity d(y, |zy—w|, 0) ts positive. 
Theorem 2. Let E satisfy the condition (H). Then for every point pyC EX 
the limit 
(12) lim , P(p; 4) = exp[Af (Po) | 


; b> bo C Ej 
exists. 

Proof. 1. Since /(p) is continuous, there exists a sphere K(y, 7) such that 
1(p) > (Po) —e, €>0, for PEE CK (ho, 7). Since po € Ej, there exists an extremal 
point q=t" which belongs to the set E K(f), 7). Therefore 


2 oi A) 2 TT | 2% exp {aL (bs) — I}. 


kj 


and 
2 D! (bo; 9 (A))| = M (ho; 7, 9 (A) exp {mA [F(ho) — e]}, 
j= 


where M(9; 7, 9 (A)) is the maximum of the absolute values of those polynomials 
in the collection 


n 


P Ik 


; (pO eee ey 
qj Tk 


k=0 
hej 


for which the point q; belongs to EX K(p, 7). 
Smices 


lim '//M (fy; 7, g(a) 21—o(7), o(v)—>0O when 7>0, 


we have 

(13) P(po; A) Sexp {Alf (bo) — el} (1-2 (7)). 
The function ®(p; A) is lower semicontinuous; hence 

(14) lim | B(p; 2) Zexp {Af (ho)}- 


1 See [17], where the proof in the case of one variable is given. 
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From (13) and (11) we obtain 
(15) P(py; 2) =exp {Al (bo)}, PoC Ez 

2. Let K(py,7) be a sphere such that /(p)</(po)+7 for PEEOK (fy, 7). 
From (7) results the inequality 


|P (6; (A) Sexp {m ATH (bo) +n]}, — PEEMK (2057). 


The absolute values of the polynomials 


BD! (p; g™ (A)) exp {— nA[f (bo) +n]} 


are therefore <1 on EOK(fy,7). From the condition (H) and from the theorem 
of C. LosTER (see [18]) it follows that there exists a number N and a neighbor- 
hood A of the point p) such that for n= N and pcA 


[2 (pb; (A))| exp{— Alf (bo) tu}}< (1 +2)", 1=0,1,-...0- 
Therefore 


2 | PY (65g (A))| SU +m) (1+ e)"exp{nALf (bo) +n]}, PEA. 
From this last inequality it follows (see [70]) that 


(16) Tim ®(p; 2) < exp {A/ (b0)}- 
(12) follows from (14) and (16). Q.E.D. 

Remarks. 1. We denote by S, the set of points cpy, where fp is an arbitrary 
point of E andc an arbitrary complex number. The function ®(; A) is continuous 
in S, if it is continuous on £, because (see Property 3) 


slim, B(; 4) = lim P(cp; A) =|c | im D(p; A) =|c| D(po; A) = O(c hy; A). 


Since C,=S,+S, it follows that ®(f; A) is continuous at every finite point p 
when E satisfies the condition (H) and exp /(p) =| Q,( p)|r. In fact, ®(p; A) is 
equal to exp {At (p)} on E (see Property 7), and cl the lower semicontinuity of 


P(p; A) it foll h 
(p; A) it follows that lim @(p; 2) = exp {Af ( (po)}- 
p—>poCcE 


On the other hand since E satisfies the condition (H), 


lim @(p; A) <exp {Af ( (Po)} 


p> po CE 


(see the second part of the proof of Theorem 2). 

2. Let E be a subset of E which contains all points of E such that none of the 
points ch, PEE, c>1 belong to E. We shall call E the boundary of the set E 
with respect to the origin (0, 0). We have (see [16]) 

d(E,|zn —Cw|,0) =a(E, 24 —C w|, 0), 


and the function ®(f; 0) constructed for the set E is identical with the cor- 
responding function constructed for FE. 
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If £ is contained in the set! {|Q,(p)|—=c>0} where Q,(f) is an arbitrary 
homogeneous polynomial of degree k, and if E satisfies the condition (H), then 
(pf; 0) is continuous at every finite point and ®(p; 0)=1 for PCE. 


Theorem 3. For every finite point p+0 the limit 


@(o; 4) ]2 


(17) Decoy? =P): 


exists. 
Proof, Let 0<A,<A, and 


(")) P Pr ze 
LY (p; p =I ae HO, AGE argh, 


re 


We denote by p” (A,) a system of points of E such that 
e+ ®,(b; 4) > max | Py) (b; Pe(”r))| 
for an arbitrary number ¢>0. On the other hand we have 
P,(B; Ag) Smax| OPP; g(Aa))| =| Py (6 Gs))| =[Ln (0; @(A2)) | exp [ Aa f (g.)], 


where g™) (A) denotes the extremal system of points. 
Using the interpolation formula of Lagrange we obtain 


®,(p a) SD) LY? (b; PP (A,))| | LO (Bj; a (As)) | exp [7 As f(q,)] 


a Ast a Aa 
<Dexp[nds/ (p Pi) ILMB; BEA) =D | OPH; HOA) LPB; LAN)! ; 
(18) 7=0 9 j= 
S(m +1) max] PH; bP) | max [LI ; pe") 

] 7 


me 


da = Ay 
<[®q(B; A) +e] (m +4) max [LY(; BEAN 


Let g™(0) be the n™ extremal system of points corresponding to A=0; then 


®,(p; 0) Smax [LNG g0))| =|L(65 10))| SDI Ew; WO)|[La( 5 be ad) 
(n+1 )max|L(p; b2(A))] 

qj 
and therefore 


(19) @(p; 0) S lim n max |Ly) (p; p™(A,))] . 


n— CO 


IA 


Since 1/2 is <0, it follows from (18) and (19) that 


1 


Ag d _ 4a 
D(p; dg) <[P(p; Ay): [P(p70))* 
Therefore 
Dp; Ay) (p; hy) 
ao) Ae Ae 


1 When Ec{|Q;(p)| =c>0} then E=E. 
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From (20) and from the Property 9 follows the existence of the limit (17). 
The function Y/(p) satisfies the inequalities 


expm <¥(p) SexpM, p= 0: 
Suppose £ is contained in the set {|Q; (p)| =c> o} and satisfies the condition 
(H). The function 


1 
2M 2 


D(p; 0) 


ea 


is continuous at every point +0. At every point py € Ej we have P(Po; =A 
and D(p,; 4) =exp[Af (po) |. Therefore (po) =exp/ (po). This equality holds at 
every point ~)€ >| E*. In fact, suppose py € £7; then 

A>0 


+ 


P(bo) = [Gite [a = exp (pr). 


. i 
On the other hand ee | <exp/(f,) for all A>0O and therefore 
0? 


P (po) S exp f (ho): 


On every analytic plane which passes through the origin and an arbitrary 
point ~, € >; EX we have 
A>0 


P(p) =exp f (by) (=const.), p+. 


Remark. In the generalized Nevanlinna theory, as developed by BERGMAN, 
the main problem which arises consists in defining functions of extended class 
(see the Introduction) and investigating their properties. 

In the present paper we showed that the classical results obtained previously 
by Lrja permit a new construction of the extended class based on the method 
of extremal points. 

One can show that use of this method provides new possibilities for deriving 
bounds for functions of BERGMAN’s extended class. In this way our results can 
be used to obtain bounds for meromorphic functions of two complex variables. 
One obtains relations between the growth of analytic functions of two variables 
and properties of their a-surfaces (7.e., surfaces where the analytic function 
f(z, w) assumes the value f(z, w) =a). 


This paper was done under a contract between the U.S.National Science Foun- 
dation and Stanford University. 
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Zur numerischen Behandlung von Anfangswertproblemen 
partieller hyperbolischer Differentialgleichungen 
zweiter Ordnung in zwei unabhdngigen Veranderlichen 
I. Das charakteristische Anfangswertproblem™ 


W. TORNIG 


Vorgelegt von W. TOLLMIEN 


1. Einleitung 


Fiir die numerische Behandlung von Anfangswertproblemen gewdhnlicher 
Differentialgleichungen sind die sogenannten Differenzenschemaverfahren bekannt 
und bewdhrt. Sie sind unter den gegenwartig bekannten Verfahren wohl im 
allgemeinen die genauesten. Der einfache Grundgedanke aller Differenzen- 
schemaverfahren ist die Umwandlung der Differentialgleichung mit den Anfangs- 
bedingungen in Integralformen. Die Integranden dieser Integralformen werden 
dann durch Polynome approximiert [4]. 

Der Gedanke liegt nahe, diese Methode auf Anfangswertprobleme partieller 
Differentialgleichungen zu tibertragen, wobei man in bestimmten Fallen und 
unter entsprechenden Voraussetzungen erwarten darf, daB hierdurch ahnlich gute 
Ergebnisse erzielt werden kénnen. Im folgenden wird dieser Gedanke auf Anfangs- 
wertprobleme der Differentialgleichung 

Oz 


Ox Oy 


= t(x, Md Oe 7 
Lu ORL OS) 

iibertragen. Eine Darstellung in der vorliegenden Form scheint bisher noch 

nicht vorhanden zu sein. In diesem I. Teil der Arbeit soll das charakteristische 

Anfangswertproblem behandelt werden, wahrend die numerische Behandlung des 

Cauchy-Problems dem II. Teil vorbehalten bleiben soll. Alle Betrachtungen 

werden ausschlieBlich im Reellen durchgefithrt. 

Vor einiger Zeit hat J. Diaz [9] das charakteristische Anfangswertproblem 
der Differentialgleichung z,,=/(%, y, 2, 2,, %,) mit Hilfe einer dem Euler-Cauchy- 
schen Polygonzugverfahren analogen Methode behandelt, wobei er auf deren 
Wert fiir die numerische Anwendung hinweist. Wahrend wir die Anfangswerte 
auf den Charakteristiken (und im II. Teil der Arbeit beim Cauchy-Problem auf 
einer Kurve) im allgemeinen in analytischer Form vorgeben miissen, ist das 
Verfahren von D1Az immer anwendbar, wenn nur diskrete Anfangswerte vorge- 
geben sind. Vom numerischen Standpunkt aus gesehen scheinen jedoch zwischen 

* 4, Teil einer gekiirzten Fassung der 1958 von der Fakultat fiir Natur- und 


Geisteswissenschaften der Bergakademie Clausthal angenommenen Dissertation des 
Verfassers; Referenten: Prof. Dr. H. K6nic und Prof. Dr. H. MENZEL. 
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dem Verfahren von Diaz und den hier mitgeteilten Verfahren vor allem hinsicht- 
lich der Genauigkeit ahnliche Unterschiede zu bestehen wie zwischen dem Euler- 
Cauchy-Verfahren und den Differenzenschemaverfahren fiir gew6hnliche Diffe- 
rentialgleichungen. 

In Anlehnung an die Differenzenschemaverfahren werden Extrapolations- 
und Interpolationsverfahren beschrieben, von denen sich die Interpolations- 
verfahren durch hdhere Genauigkeit auszeichnen. Man erhalt stets diskrete 
Lésungswerte. In den Nummern 2. bis 5. werden die Vorschriften fiir die prakti- 
sche Rechnung — Extrapolationsverfahren, Interpolationsverfahren fiir das dort 
naher beschriebene ,,Anfangsfeld“ und fiir die ,,fortlaufende Rechnung‘‘ — ent- 
wickelt. Bei den Interpolationsverfahren werden die Naherungen iterativ er- 
mittelt. Die Nummern 6. bis 8. enthalten Kriterien fiir die Konvergenz dieser 
Iterationen und das Prinzip einer Fehlerabschatzung fiir die Werte des Anfangs- 
feldes. In 9. findet sich ein Beispiel. 

Die wichtige Frage nach der numerischen Stabilitat der mitgeteilten Verfahren 
wird hier nicht beantwortet. Die Betrachtungen dariiber sollen einer spateren 
Arbeit vorbehalten bleiben. 


2. Ausgangsgleichungen 

Es sei das charakteristische Anfangswertproblem 

Oz 02 
#2, 


=e (x, V, 4, 
(1) 
2(x,n)=o(x), 2(6¥)=T(y), mit o(€) =7(n) 
zur Losung vorgelegt. Bekanntlich lassen sich alle hyperbolischen Differential- 
gleichungen 2. Ordnung in zwei unabhangigen Veranderlichen mit linearem Haupt- 
teil aut die Normalform (1) transformieren. Die praktische Durchfiihrung dieser 
Transformation wird manchmal nicht leicht sein. Darin liegt die einzige Schwie- 
rigkeit bei der Benutzung der im folgenden zu entwickelnden Naherungsverfahren. 

Die Frage nach der Existenz und Eindeutigkeit der Lésung des Problems (1) 
beantwortet der [2] [3] 


Satz 1. In dem Rechteck 
Mt awd, q=<v<d,, Z2,u,0 dehebig, 


sei die Funktion f(x, v, 2, u,v) stetig, in jedem abgeschlossenen Tetlgebiet g von K 
beschrankt und erfiille dort eine Lipschitz-Bedingung 


| f(x, V7, 29, Ue, Ve) — f(*, ¥, Z, Uy, V4)| = M(|z,— 4| a | %. — u| a Woes v,|). 


Es seien in den Intervallen a,<x<b,, <y<d, ferner die stetig differenzierbaren 
Funktionen o(x) und t(y) vorgegeben und (&,) ser ein innerer Punkt von KR mit 
oa (£)=t1(n). Dann hat das Problem (1) in K die eindeutig bestimmte Lésung z=y (x, y) 
mit p(x,)=a(x) und y(&, y)=T(y). 

Es ist zu beachten, daB durch die Vorgabe der Anfangswerte o(x) und t(y) 
auf den Charakteristiken y=7 und x=& die Lésung des Problems schon ein- 
deutig bestimmt ist. Man erhalt sofort 

op(*,n) _ oy(s,y) _ oy 
WM — g(x), WEP) — ay). 
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Die Differentialgleichung (1) stellt lings *«=é eine gewoéhnliche Differential- 
gleichung fiir ee! dar, da only dort bekannt ist. Ebenso erhalt man eine 
x 

gewohnliche Differentialgleichung fiir —_ langs y=7. Die Losungen dieser 
Gleichungen existieren und sind eindeutig. Damit sind die z, z,, 2, auf den An- 
fangsgeraden festgelegt. 

Das Problem (1) geniige den Voraussetzungen des Satz 1. Mit den iiblichen 
Abkiirzungen z,,=s, 2,=, 2,=g ergibt sich als Lésung: 


z(x,y) =2(x,n) +2(E, 9) — 2(E7) +f it v,2, p,q) dvdu, 
p(x, 9) =pl(x.n) +f fle 0,2, 6.9) a0, (2) 


q(x, 9) =9(§,») +f tu, Y, 2, p,q) du. 


Die Integranden in den letzten beiden Gleichungen sind Veranderliche nur von 
x allein bzw. von y allein. 
Die Intervalle a4,<aSxSb<b,, cy<cSySd<d, werden unterteilt durch 
A= XH XyQ<--- << 4K; =), 3) 
C= V9 Wy Yn yy =a. 
Zieht man durch die Teilpunkte der Intervalle Parallelen zu den Koordinaten- 
achsen, so wird das Rechteck al xb, cS y<d inj / Teilrechtecke unterteilt. 


3. Extrapolationsverfahren 
Die Rechnung sei schon bis zu Abszissen x=x,, y=y, fortgeschritten und 
daher Naherungen 


2 (%,,, Vy) =2u ys D (Xu Vy) = Pury G (Xu Vy) = Qu,» 
f (Xu, Vy» Suv» Pig. Cae Pa be 
(Was Opting 02 Ween On Ae es eee eecar | 
bekannt. Es wird nun in (2) f(x, y, 2,6, q) durch das Polynom ye CP Ae 
das an den Stellen (x,_,, y,_,) die Werte | ns Tommasi 
Pe OA es 0) Mt Se WES), 


Sei stets bei allen folgenden Betrachtungen 
Agree hy hs Vie po Voce (fee 0515s faethy a0, dance l eel 


dann ist mit der Transformation [7] 


. 2 UH yn : 
Ug Css y) = > ao Shih o Ve V; Pass (4a) 
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Es wird ferner {(x,,1, ¥, 2, ,@) durch das Polynom P,(y) ersetzt, das an den 
Stellen (%,;1, Ys_,) die Werte /,,,; ,_, annimmt, schlieBlich f(x, y,.1, 2, p,q) durch 


das Polynom P(x), das an den Stellen (%,_,,, Y,41) die Werte /,_,, .4, annimmt: 


— Be ue ee 
yA 5s ts —Vyf Var S$? (4b) 
= ide ae Ht) 
Be Gear: Ve if Sek" (4c) 
Nach (2) erhalt man dann mit ae "=, die Naherungswerte: 
By Teese == Ay Ty, 6 IE |p, SAA, Zz, .thk m2 oi BB, Ve Ue 
p=0) v=0 
Prorsti= Pris th 2B, Ve Teas) 
(5) 
Vy+4, Sera ==". Sed +h > Be Ve ibs Soe 
u=0 
Pee sti=H(%41. Vsta> Bere pas Praesens Geass 1) 
mit 1 
ag fred w= f u(u tt) (+2)... («+e —1) da. (6) 


0 


Tabelle 1 [5]. Gewchte B, 


oom 0 | 1 2 | 3 | 4 5 
P | 1440 720 600. 540 025 475 
e | 4440 1440 | 1440 | 4440 1440 1440 


Setzt man /(x, y, 2, p, 7) =F (x, y), so gilt fiir die exakte Lésung der Differential- 


gleichung: 
2(%,44 ) Vs41) J 2(%,41, ys) a z(x,, ¥541) a z(%,, ys) 2 


ar hk De > B.By Vi VEX, ys) f. Sdn eL 
p=0 v=0 (7) 


P(%p44> Vo+1) = P(%pa1, Vs) att k XB, Vy PX a1, ys) ab Sua 


m 


q(x Pag ihe Vou) = ee q(%Xp, Vs41) a8 h a lay Ve Ta Vo41) a Spe i 
u=0 
Die Kubatur- und Quadraturfehler kénnen abgeschatzt werden [J]: 
pp entt 
Alga +hk si Been “Oyn tt Fe 
Omrny2 


oxmti ayntt 


om. 


+ 


rene Rs Pie Bacal Oxmt1 


as er ater |Bu+2 


eee 
onrl (8) 
Sane =< Re nelle imc oynti ee 
m+1 
| Sodesabe Bast a Re 
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Mit den Formeln (5) lassen sich fortlaufend neue Naherungswerte der Lésung 
berechnen. Da die praktische Differenzenbildung nach 2 Verdnderlichen recht 
umstindlich ist, erscheint es zweckmaBig, statt der Differenzen direkt die Funk- 
tionswerte zu schreiben. Aus (5) erhalt man dann nach kurzer Rechnung: 


m 
4r+1,st1 = 441, By 542 a a a hk Z 
p= 


n 
Pyar eat = py +1,s i k Ls n,2br+t,s—v> 
v=0, 


m 
Gy+1,s+1 — Uy, s+1 +h Xm brn se 
iA 
beaten Ii 7 as Pee Oe 
H : A 
a= (— 1) > B,(7)- (10) 
Pay: 


Fiir die exakte Lésung der Differentialgleichung gelten entsprechende Glei- 
chungen mit den Abschatzungen (8) fiir die Kubatur- und Quadraturfehler. 


mit den Gewichten 


Tabelle 2 [5]. Gewtchte a; 1 


i] 


0 1 radia fe Ge ted i ee 
ee 4277 
2 12 24 720 1440 
, 1 16 Zdieeam Wn bags re 6 
2 {hal - 24uc) “98 1440 
x Srey her By 2616 9982 
12 24 | 720 1440 
; hin 9.9 heehee 4 7208 
24 | “726 1440 
" 251 2877 
| 720 1440 
5 SeieN, 
1440 


Die Approximation der Integranden in (2) kann offenbar auf sehr verschiedene 
Arten vorgenommen werden. Es soll darauf aber nicht naher eingegangen werden, 
zumal durch Extrapolation im allgemeinen nur ,,Rohwerte“ gefunden werden 
sollen, die dann durch Interpolation verbessert werden. 


4. Berechnung des Anfangsfeldes 


Um die Formeln aus 3. und die im nachsten Abschnitt 5. entwickelten Inter- 
polationsformeln anwenden zu kénnen, benétigt man schon entlang der Anfangs- 
geraden ein ,,Anfangsfeld‘‘ von Werten. Zu seiner Berechnung kénnen alle 
bekannten, hinreichend genauen Verfahren dienen, wie z.B. die verbesserten 
Differenzenverfahren und Reihenentwicklungen. Im folgenden wird ein Ver- 
fahren vorgeschlagen, das unter den Voraussetzungen des Satz1 in 2. stets 
anwendbar ist. 
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Auf den Charakteristiken x=x) und y= yp seien die Anfangswerte z(%9, ¥) 
und z(%, ¥9) vorgegeben. Dann kennt man auch die f und gq langs der Anfangs- 
geraden. Dabei sind £(%9, y) und g(x, yo) eindeutige Lésungen gewéhnlicher 
Differentialgleichungen. Gelingt es nicht, diese Gleichungen exakt zu lésen, 
dann empfehlen sich hier zur numerischen Lisung besonders die Differenzen- 
schemaverfahren. (Siehe z.B. [5]). Man verschafft sich nun zuerst Rohwerte 
2°), pl), gf), f°), indem man z.B. die Formeln (9) benutzt und 7—1 statt m, 
Bay) Salts SetZt, (ety 2h ns, MESA, 2005) 


v—1 s—1 
a 0 
Bl gS ate, ey or eR D, nee 1,u%s vad aires 1—ys 
=0 »=0 
Ss i ; 
0 0 
(De je em k 2 ee ee: 
v=0 
r—1 (41) 
0 
Ge aie +: h > Cet epee 
u=0 
ete re ees Pee 
Dabei ist 
219) = z(%;, y1) 
[0] __ p(x. : 
Pit P(x), 1) Tin) ef O godet? 20 oder 7 =1= 0: (12) 


qj,i ay q(%,, y1) 
i) = F(x,, y1) 


Die so berechneten Rohwerte werden nun durch ein Interpolationsverfahren ver- 
bessert: Mit [Z] 


2 EG 
P (2) sul (ene PPE) = nes 
ist 
bw Vea) 
Bowles 2, 2, BPm(x,) B™ (yy) lee (13) 
u=0 y= 


das Polynom, das an den Stellen (x,, y,) die Werte 7, , annimmt, (uw=0, 1, ...,™; 
y=0,1,...,”). Mit der Transformation 


H=Xpuh, yvy=V%+tvR 
wird daraus 


Lea SUSE UD (n) 
Pl 9) = Oma) = YY Gita “oo fm (14) 
mit 4 
0 (wv) == J w —») 
v=0 
Ferner ist ™ aa. 
i) = ia te (15a) 
v0)” 
das Polynom, das an den Stellen (x,, y,) die Werte /, , und 
— eae (m) 
Qn (H) = > SE has (15) 
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das Polynom, das an den Stellen («,, y,) die Werte /,,, annimmt. Setzt man in 


(2) =X), N=, ersetzt die Integranden durch die Polynome (14) und (15), 
so erhalt man die Interpolationsformeln : 


m on 
5 n 
Bp é Pr, 0 eF0,4hy Court PBs PER Re ah IW 


u=0 »v=0 
Drs =F Pro hi k a Mavi iy: 
Ge (16) 
Gs = qo, s ie h ve Vis ne 
p=0 
hes=T(%,, Vexbg rcs Preteen (7 


mit 

Cy meme iea e er oe Py, 

Vee A ee dw a | oe 1)...(w—A+1) (w—A—1) ...(w—7) dw. (417) 
0 0 


Tabelle 3 [6]. Gewichte y!, , 
Die Zahlenwerte fiir die Matrizen 


7 ft ; 
Misc Gases . 


: 
ra 
Yo0 Ve1-°° 


a Vig me Ve 
sind: 
ale 
at as 
oegem 
J=2 
8 1 
cries be ee 
fy S24 1 ‘ 
3 3 3 
j=3: 
9 19 5 1 
2h BA | Boe 
n 4 1 . 
3 3 3 
3 9 9 3 
8 8 8 ee 


Fir die exakte Lésung der Differentialgleichung gilt: 


2(%,, ys) =2(%,, Yo) +2 (Xo, Ms) == 2Xyy, Vo) +hk yz > Gyorgy ll Coe Vee iy ate 
uH=0 v=0 


P(%rs Vs) = b(%,, Yo) + kD V3, F(x,, Y,) + Shea, 
v=0 


(18) 
G(s) =9(%0, Vs) +h Dy, Flay, Ve) + S84. 
0 


Y= 
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Die Abschatzung der Kubatur- und Quadraturfehler ist hier unbequem; man 
erhalt nach einiger Rechnung [7]: 


iS m4+2 m-+-1 
[Sriantal Sskh ae afl 5 lyin m+1 Viem+1| 7 
a +1 
ale n+2| o n+1 n+1 
rhk aynti |. an lytt nea — yeti aii 
* 1=0 
1 fymr2 pre naa 


OxM+F1 Qynt+1 


y—1 s—l1 
m+1 

” + [et mti— Vy met | |Yeeinea— Vener (19) 
o =0N/=0 
a s—1 
Oy = | n+1 

my eR ail, = 
Max [0 


rove 8) are ara 


emt. 


= r—1 
Cr, s m+ 2 m+1 m+1 
Srna | Sh sore Di ed x Vim+1| 2 
x Max f= 


Es sei vorausgesetzt, daB die in (19) auftretenden Ableitungen beschrankt 
sind. Dann ist mit k=h, m=n und fiir h-0: 


Santi) = oh"), — Shis(h)= ow"), — Shs (h) = 0(h™*). 


Nun zur Anwendung der Formeln (16): Man kennt bereits Rohwerte 2!?},, pi), qi?) , 
(0] 


wy: Erste Naherungen werden dann berechnet mit Hilfe der Iterationsvor- 
schriften: 


Sy 


Piss by, 071 0,40 ee Sym at, fl), 


p=0 v=0 


pe. =Pro th DM, 


Ge Lag, bh LMA, 


u=0 
ea | Cnn ee We ee CAs Sa Oe Nn), 


mit der Schreibweise (12). 


Bei der Berechnung weiterer Naherungen vereinfachen sich die Vorschriften 
(20): Die (A+-1)-ten Naherungen sind: 


Pint se zi] +hk Sy S on Wye A (ff) — (BG 1), 
1 


ce 


n 
ast ret eta Dew beaee 
(21) 
Ge = Gna h 3 Yiullins — fins) 


(tier AC Se y, coe Wags : gh) ; 
Ftp els 
Hinreichende Bedingungen fiir die Konvergenz dieses Verfahrens werden in 
6. gegeben. 


Arch. Rational Mech. Anal., Vol. 4 29 
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Mit Hilfe der Formeln (11), (20) und (21) berechnet man zunachst Naherungs- 
werte tiber dem in Abb. 1 (fiir m=n= 2) mit ,,1‘° bezeichneten Rechteck. Um 
die Verfahren der fortlaufenden Rechnung anwenden zu kénnen, bendtigt man 


noch Werte iiber den in Abb. 1 langs y=, mit ,,2“, ,,3“,... und langs x= %9 
Y4 mit ,,2’‘, ,,3’‘, .... bezeichneten Rechtecken. 
Ya | = Die Formeln (11), (20) und (21) dienen zur 
| Berechnung auch dieser Werte, man _ hat 
lediglich 
Y6 i 5 
| bei der Rechnung langs yy, fiir die z, ,, 
; | PD 4s Garay den Index o durch 
: | Mm+o, (Mie 4 2) .6 oho == 12ers Oe 


bei der Rechnung langs x= x, fir die 2, ,, 
Poa» Io,a» fo,, den Index A durch 


a Cero 
yea eset eee oe 


Yor, e Niwa A, (IN S47 A = ease 


Abb. 4, Zur perce des i one zu ersetzen. 


5. Interpolationsverfahren fiir die fortlaufende Rechnung 
Man ersetzt hier /(x, y, z, f, g) in (2) durch das Polynom P*,,(x, y), das an 


den Stellen (%, +1.) Vs+1—»)) (U=0,1,...,;7=0,1,...,), die Werte fis DRG, 
annimmt. Entsprechend werden f(x, v, z, p, q) ae ie yreipy gy. ersetzt. Da 
die Herleitung der Interpolationsformeln ganz analog der Herleitung von (9) 


erfolgt, konnen wir uns sehr kurz fassen. Mit der Transformation 


a VAD Soe 
h k 
und =4%,, 7=y, erhalt man zunachst: 
“ 4 
ae i eS * 
By Ro ope ae Se Bo pag ss ye Vy f r+1,s+1> 
p=0 »v=0 


Pr4t,s4+1 = Pr i ee k ve Vy bce s+1> 


m (22) 
Cp seen = UP a) h >» 15 Vi Eas 
ie 1s¢1 —T(% 44) Vera 2p4 ee es oe Qy-+1, s-+1) > 
mit 1 
Aik 
fis. (u—1)u(u+1)...(u+1— 2) du. (23) 
0 
Tabelle 4 [5]. Gewichte Bt 
r | 0) | 1 | 2 | 3 | 4 | 5 
6 | 720 | 360 | _ 60 30. 19 27 
| F250. 950) aaa 720 720 720 


Die in (22) auftretenden Differenzen werden wieder durch Funktionswerte 
ersetzt. Auf der rechten Seite der dann entstehenden Vorschriften kennt man 
den Wert f,,1,,, nicht, sondern nur einen Rohwert pe 41. Man wird daher 
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die Naherungen wieder iterativ ermitteln mit Hilfe der Gleichungen: 


(1) a 
Orel ait —— Syl 6 Teer epee ey, 6 


m n 

+ io Sans ie TS 2 a CO des n, aa BL, oa) 
u=0 »=0 
Ait 


pu, s+1 —_ = Py 41 s = k lt Xp, pee s+1 a > aye v pan ant ? (24) 


y=1 


1 * 4[0 x 
Greta 97 Sei h ee afi Ste p mw penal ) 


Par oes hee Ned) Vs41.# oon ras tagglani)s 


mit 
= (1 (? )Br. (ZO) 1 e2sR ey, (25) 
A= 
Fiir A=1 vereinfachen sich diese Formeln betrachtlich: 
y= ear WR 0 no fae — a) 
Wai Ses oes SE sf k Xn, poe, ea hase ea) ? (26) 
qQtaa SE Gena + hoo (fy, 44 — Wert cea); 
ees =f (x, t1» Vs+1> a 11» ey > qe eile 
Tabelle 5 [5]. Gewichte aj"; 
Fiir die exakte Lésung der Differentialgleichung gilt: 
2(%,41, Vo41) = z(%, 41, ys) i z(%,, are) = 2(%,, ys) ar 
m n 
i h k > o3 CCE OCA Vora oi) oe Charrette 
u=0 »v=0 _ (27) 
P (X41 Vs41) =P (% p41 Vs) +A M2 Op eT (X44; 6 eet ae | Saas 
yv=0 


q (Xp415 Ys41) = q(%,, Vs41) -F h ¥ cos, 7 TG ea Veit) + Sia: 
2 


29% 
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Abschatzungen fiir die Kubatur- und Quadraturfehler sind [7]: 


epee en = k jae ee +1 


+1 
ie ans set Be A 


oxmth Saar 
m+2 4n+2 Geta s 
“tr h k |Bn+1| |Br+a| Oxm+i FL Gynt Mane’ 
ax (28) 

onrt 
Se Piss Wee oes 41 Gynt Maxi 
— em rt 
[Shea] SH" [Baal | ees P| 


Den in (24) bei der ersten Iteration bendtigten Wert //°/,.,1 kann man sich mit 
Hilfe der Extrapolationsformeln (9) verschaffen, oft geniigt aber auch eine 
Schitzung dieses Wertes. Man wird auch wegen der geringen Rechenarbeit 
einige Iterationen mehr in Kauf nehmen, anstatt vorher durch die etwas um- 
standliche Extrapolation einen besseren Ausgangswert fe, ..1 Zu berechnen. 

Dem ,,Verfahren der zentralen Differenzen“ fiir gew6hnliche Differential- 
gleichungen [5] entspricht folgendes Interpolationsverfahren, das wiederum als 
Iterationsvorschrift benutzt wird: 


(1) ee y 
Spe, SLU Spal etary epee at cee 


+S stat fasta tb haasa t fete + 
=e A Ve 44 iyo Tges fea his) Ol ee 


1] R_ +410] ee 
P41, s+1— = br41,5-1 -- cy Re t1,s+1 + he t1,s—1 = Af, 5] , 
1 h -sto 
| ie 1=—% Lsaek a rh Peres os ee con =F 47,3] ’ 
1 f: 0 0 
fina > (Kya, YVs+1> ay Bi , pies Gn Sa) . 
Fiir A=1 vereinfachen sich die Gleichungen zu 
Fatal Ee A(t | A 
ata = Area, s41 + ( OD Se Bs Ay 
acy wai k i ata = 
b; aie b pb} A s+1 + — an ralizonrs! aa paaty) ) 
LA [4] [A—1] (30) 
(bi OM nits == se Reese ee faa. his 
[A+1] [A+] [A+1] 


A+1 ee 
Tani ead = f(x. Patctlay) Vsti» % y+1,s+1> pan s+l> G74 1,a4s1) . 


Fiir die exakte Losung der Differentialgleichung gelten entsprechende Gleichungen 


(vgl. (27)) mit dem Kubaturfehler S¥¥ und den Quadraturfehlern S** und S¥*, 
die wie folgt abgeschatzt werden kénnen: 


[she ieee o4 a5 kh Cn hei? oh 
y Ab Owe Max 45 | oy* Max 8100 | dx* dy4 Max ” 
fe 1 
Cake AE seep < 04 (31) 
90 oy* Max ‘ ~ 90 | x4 Max ; 
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6. Konvergenz der Anfangsiteration 


Beriicksichtigt man, da8 die Differentialgleichung (1) eine Lipschitz-Bedingung 
mit der Lipschitz-Konstanten M erfiillt und setzt man 


tomer Be ET Me PE P= EL, aM — A= BE, 2) 
[5 | + [OS | + 1802 | = 08, 
so gelten fiir die Gleichungen (21) die Abschatzungen 
ja sham SD yal bybel wha 
| 64| < <pbM x ly, | wri, 


[5a] saat S toe | wi 
Gel 2 eat aS A Os) 


Die Addition dieser Gleichungen liefert 


wi<hkMy Div ol [ool we +4 D [yt niger es M4 hM Sot wlts a (34) 


u=1 v=1 p=1 


Dieses Ungleichungssystem hat die Form 


ws> Says] ws. (35) 


eee a 


X= ([4z31) 


dem Betrage nach <1, so ist die Konvergenz der Reihen 


Sind alle Eigenwerte der Matrix 


oe) A CON 
D coe Dioley S64, Abel, mis Si, aha); 


A=0 A=0 A=0 


gesichert. Daraus ergibt sich die Forderung 


Max { > x >» [Ar Jh<s (Zeilensummenkriterium) , 
oder te (36) 
Max { > 2 Dy [Ag sf <1 (Spaltensummenkriterium) . 


Hieraus erhalt man obere Schranken fiir die Schrittweiten # und k. Fir k=h 
m—=n=2 lauten diese Bedingungen 


=| 720 ) 
ee tS eap aie we 


a it —1). 


bzw. 
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7. Konvergenz der Iterationen fiir die fortlaufende Rechnung 
Mit den Bezeichnungen (32) gilt fiir die Gleichungen (26) die Abschatzung: 
wil, oe < (hk| op, o| | oe, o| = R| on o| + h| on, o|) wiry EY: 
Die Konvergenz der Reihen 
Di OF aa py él, Sip si da, s+1 
4=0 
ist daher gesichert, wenn 
M(h k| ach, o| |o%7,0| + & | %n,0] + 1 |%m,0]) <4 (38) 
Mit k=/A und m=m erhalt man: 


h< a (|! bo 1). 


ZahlenmaBig lautet diese Bedingung: 


n<2(/1+4—3) fir m=1, 


M 
h< = (|/ lepers 1| UA — 
Wh me (= 7 fir m=3, (39) 
h< os Vi ae 7 Uae 
yee 88 (n+ 14 11] fir m=5. 
Fiir das Verfahren (30) toe die Rae ace 
pos(r 4) ao) 


8. Fehlerabschatzung fiir die Werte des Anfangsfeldes 

Es werde vorausgesetzt, daB die Iterationen (21) bis zum Stillstand durch- 
gefiihrt wurden, so daB die Indizes [A] fortgelassen werden kénnen. Die Ab- 
rundungsfehler sollen auBer acht gelassen werden. 

Es wird angenommen, da die Berechnung des Anfangsfeldes entlang der 
Geraden y= yy bis zu einer Stelle x= x, fortgeschritten ist. Es sind daher bereits 
Werte 2,45 Du iuiQarve fan (MeO, donne es PaO, 4s op ey, bekannte Denktanan 
sich die Anfangswerte auf den Geraden x=x,, v=Yg vorgegeben, so lauten die 
Formeln (16) zur Berechnung von Naherungswerten 2,,,5, Petys> Iatrso f 
(2=0, 71, 2%, ....)7 


a+r,S? 


m n 


= | | mm n 
Zatr,s = Sat+r,0 | Za, s 24,0 | hk > pay More ati se ao 


” R=0=0 
Powers on Date a k 2 is, v Lines n? 


Vaty,s = = Ga, s oF h = Print a+b, s? 


u=0 


ers Cet) Vs» Zatnss Derrlsn Tee)? aed, Phy cong hes Seay 2 ry It) 
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Dabei sind alle Werte auf der Geraden y=y, exakt vorgegeben, wahrend die 
Werte auf «=, bereits durch dasselbe Verfahren berechnet wurden und daher 
mit einem Fehler behaftet sind. 


Fiir die exakte Lésung der Differentialgleichung gilt: 
B(Xg+y, ys) a B(Xatys Yo) ‘ir R(X, ys) i R (Xe, Yo) at 


=i hk a ye NEO ca sn (Sianatt yy) 7 SET n+1> 
u—Ow—0 (42) 


P Gacy ys) == Pa 453 Vo) aie k yy Vey l fay Vy) =f HAG 
v=0 


Y(Xuty A) are I (Xa, 4) = (Ps Mee. Nasty’ ;) a Ssh. 


u=0 
Setzt man fiir / =>1 


Rj — 2(%;,M) =F Pj1— P(%;, 1) = & 1, 9,1 — V(%j. V1) =e 


subtrahiert die Gleichungen (42) von den entsprechenden Gleichungen (41), so 
erhalt man mit der Lipschitz-Bedingung die Abschatzungen: 


, 


| exe )5| OED, 2 lyre lyz,| {| Epiu ol +|Exsu,»| ae Veal S Rie. 
Bal v= 


n 
Les yas| ae ee {| ae - \Eeasr alate le ae. [Rize s 
= 


(43) 


Bese AED [rel {Sesmel +Barmel + etn lS Rese 
mit . 
Ratr sl = [0,0 +e-AMD |yttol [rel {eel +1821 + Eee b+ [Stats 
[Ra+r,s|=[Srial, (44) 
[Retr sl= [8,5] +4 M|y%ol {lex 6] +18,1-+ 18, ]}+1 Sea], 


Gree A 2s NaS ale em: 


’ 


Die rechte Seite dieses linearen inhomogenen Ungleichungssystems ist bekannt. 
Die Matrix des Systems, hier kurz mit D=(d;,) bezeichnet, ist von ,,monotoner 
Art‘ [5], denn das hinreichende Kriterium hierfiir ist erfiillt: 

1. Es ist d,;, SO fiir ¢=ER. 

2. D zerfallt nicht. 

3. Es existieren t>>0 und 3 >0 mit 23=t. 

Die 1. Bedingung ist immer erfiillbar bei geeigneter Reihenfolge der Summa- 
tion, 2. ist stets erfiillt. Um zu zeigen, daB auch die 3. Bedingung erfiillt ist, 
werden die Ungleichungen (43) majorisiert: 


Vay We +a as 2 lyre 7%. ares Ve ae PY an a Tey 
p=1 v= 


Zope rs kM >, Bae (Aner ae po Dyas) a Dre, 
v=1 


Luins hM >, lana (Za aZ, tee 24.25) =F Lares 
1 


AO, BI SSA 2s he), 
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mit 
kee eal Zee Eiaree| ee 
Leesan = Laue ) Bs eates s| Ss i eee Sea he (46) 
|For, Lae bade a a+r, s 


(af a l, 2a ee WP, SE Ae 2: Sas BY: 
Das Gleichungssystem (45) hat mit © als Einheitsmatrix die Form 


p= (C Dg et (47) 
und mit E—D=°A die Lésung 
p= (ELU+ WU +---)t. (48) 


Diese Matrizenreihe konvergiert, wenn alle Eigenwerte der Matrix % dem Betrag 
nach <1 sind. Nach (36) ist das stets der Fall. Damit haben wir eine Darstellung 
der Form 

4—= Gt. 


Die Matrix & hat nur positive Elemente. Wegen t>0 (d.h. jedes Element 
von t ist >0) ist auch 3>0 und damit auch die 3. Bedingung erfiillt. 


Wegen der Monotonie des Systems (45) sind dann die Z,, ;, a Z,,., obere 
Schranken der | ¢, |, | €.,|, | &.,al- 


Sei a= Max (WU) der Betrag des groéBten Elementes der Matrix YU, so gilt fiir 
den Betrag des groBten Elementes Max (2?) der Matrix QU?: 


Max (2?) < (Max eS ax) Max a;,—=Aa, 
4 k 4, 


allgemein 
Max (l") Sa 2a 


Nach (48) erhalt man so fiir A<1 die Fehlerabschatzung: 


= 


i Q 
oa! 
Me 


Vem = T41,s +—— (Te le eee 


pel el 
= im nm wn = 
Beare alt acest ht Al 2, 2, Ferny t Tame Tesi oN (49) 
BM=1 v= 


Casas oh Tepe av ee) : 


Ms 
MR: 


R 
a 
© 
-—— 


=| 
Ifa 


= 
I 
ran 
I 
ran 


Fiir «=0 treten in diesen Ungleichungen anstelle der oe T, ,, T,, , die oberen 


Schranken (19) fiir die reinen Kubatur- und Oundeariricnic: 

Sei jetzt die Berechnung des Anfangsfeldes entlang der Geraden x=, bis 
zur Stelle y=y, fortgeschritten, (B=n,2n,3n,...). Man erhalt dann ganz 
ahnlich fiir die Fehlerbetrage der Naherungen z, p+s» Prpt+s> 4,p+s Abschatzungen, 
die nach Betrachtung der Abschatzungen (49) gleich hingeschrieben werden 
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k6nnen: mon a 
lSors\= topes a > > es I re Te 
? oe eS 
[e,e+s|SUrpes+ “s » Chee ‘bute ils B49) > (50) 
re eet ag Be 
opel Wee, ene ltt st tapes) 
sit i sept! 


Trove = |€y,al +h k M| yo] 2 |yeul {| Gl + l,l + 18,0134] Sram sal 
oa 


Traty = |E,,0| +2M ye ol {eel +821 +1&, el} t+ |Setal, (51) 
(hep pe Keo ) 
(P=, 205 30,..2; Bee 1,2... ; S12, ee) 
Nach (45) haben die Werte a und A fiir k=/ stets die Form: 
a=bMh, 
A=BMh 
mit konstanten 6 und B, so daB 
Gigs) ae init. UB pee Mists (52) 


1-A 1—BMh 
Fiir alle Verfahren der fortlaufenden Rechnung ergeben sich wegen der Tat- 
sache, daB die Fehler nur dem Betrag nach abgeschatzt werden k6énnen, sehr 
ungiinstige Fehlerschranken, die nach wenigen Rechenschritten unbrauchbar 
werden. 
9. Beispiel 

Es wurde bewuBt eine nichtlineare Differentialgleichung ,,konstruiert“, um 

die exakte Lésung mit den Werten der Naherungslésung vergleichen zu kénnen. 


Differentialgleichung: 
Oz oe oe al 


Oety hoe Ooi e 
oder: 
ce Bus “ Simon 
Anfangsbedingungen: 
(x, 0)", 2(0;¥) = cos y. 


Daraus ergibt sich: 
p(x,0)=2x%e", 4(0,¥)=—siny, 


£(0,y)=ycosy, 9(x,0) = xe". 
Die beiden letzten Bedingungen werden bestimmt aus 


'(0,y)=—P(0,y)tgy+cosy, (0,0) =0, 
q'(%, 0) = 2q(x, 0) * +e", q(0, 0) =0. 
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PLO) Ant Se en CQ 


ice) 


0 0) 0 ©) @ 0 
w= C= bb a 6A OG 
(er vE+ = EG — oe Lo 
l= Va 9b+ Og 9L+ UC 
$9 — eis 0) Sis oi Ch 
baa em Cla OGe Mey CGR 
co (ES Si p= (= — ol 
Lae Shes b= eS Va ho 
= Or+ lEe— Ve= 6¢b— iv 
8 Z 9 S 4 e 


Qa 1s ee SSS S 


sOb- "9 Aaqyay 


| 
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Loésung des Problems: 


-2 - 
Det COSA, 


Ergebnis der Rechnung: 


Das Anfangsfeld (umrandet) wurde nach (11) und (20) bzw. (21) mit m=n=2, 


h=k=0, 1 berechnet. Fiir die fortlaufende Rechnung wurde Verfahren (29) bzw. 
(30) verwendet. Bei der Berechnung des Anfangsfeldes waren durchschnittlich 
drei, bei der fortlaufenden Rechnung nur ein bis zwei Iterationen erforderlich. 
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Zur numerischen Behandlung von Anfangswertproblemen 
partieller hyperbolischer Differentialgleichungen 
zweiter Ordnung in zwei unabhangigen Verdnderlichen 
II. Das Canchy-Problem* 


W. TORNIG 


Vorgelegt von W. TOLLMIEN 


1. Einleitung 


In diesem II. Teil der Arbeit werden numerische Verfahren zur Behandlung 
des Cauchy-Problems der Differentialgleichung 


Oz =i (x re Le sae 
BEo yp Pea Ba 


mitgeteilt. Die Anfangswerte sind dabei auf einer streng monotonen Kurve vor- 
gegeben. Nach einer passenden Koordinatentransformation kann ahnlich wie 
beim charakteristischen Anfangswertproblem vorgegangen werden. {Siehe Teil I. | 
In den Nummern 2. bis 5. werden Extrapolations- und Interpolationsverfahren 
zur Berechnung des Anfangsfeldes und fiir die fortlaufende Rechnung beschrieben 
und in 6. bis 9. theoretisch untersucht. Es gelingt, im Prinzip Fehlerabschatzun- 
gen fiir die Interpolationsverfahren anzugeben. In 10. wird fiir das Interpola- 
tionsverfahren der fortlaufenden Rechnung die Konvergenz der Naherungswerte 
gegen die Werte der exakten Losung bei kleiner werdender Schrittweite gezeigt. 
In 11. wird ein Beispiel durchgerechnet. Die Frage nach der numerischen Stabilitat 
der Verfahren wird an dieser Stelle nicht erdrtert. (Siehe Teil I der Arbeit.) Der 
Teil I dieser Arbeit wird als bekannt vorausgesetzt, die Bezeichnungen werden 
teilweise tibernommen. Wir kénnen uns daher etwas kiirzer fassen. 


2. Ausgangsgleichungen 
Vorgelegt sei das Anfangswertproblem (Cauchy-Problem) : 


s=f(x,¥,2,0,49), 
a(x, p(x))=a(x), p(x, p(x))=B(x), a(x, g(%)) =y (x) 


auf der Anfangskurve y= @ (x). 


(1) 


ie BPs Teil einer gektirzten Fassung der 1958 von der Fakultat fiir Natur- und 
Geisteswissenschaften der Bergakademie Clausthal angenommenen Dissertation des 
Verfassers; Referenten: Prof. Dr. H. Konic und Prof. Dr. H. MEnzeEt. 
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Die Frage nach der Existenz und Eindeutigkeit der Lésung dieses Problems 
beantwortet der [2], [3] 
Satz 1. In dem Rechteck 


Ms By So Op, _ 2) 4, 0 beltedig, 


set die Funktion f(x, y, 2, u, v) stetig, in jedem abgeschlossenen Teilgebiet g von KR 
beschrankt und erfiille dort eine Lipschitz-Bedingung: 


[7 (%, Y, 22, Ma, V2) — f(%, Ys 21, My, %4)| SM (|2q — 4| + |g — m%| +] ¥2— 241). 
TES Set 

RK: y= p(x) mt 1<x< 4%, eine Kurve, die ganz in K verlauft und die stetig 
differenzierbar und streng monoton ist. Dann hat das Problem (1) in KR das ein- 
deutig bestimmte Integral z= (x, y) mit p(x, p (x)) =a (x). 


Die # und g konnen auf der Anfangskurve beliebig vorgegeben werden, miissen 
jedoch der Streifenbedingung 


dz—pdx+qdy (2) 
gentugen. 


Das Problem (1) wird modifiziert, indem die neuen Koordinaten 
ce pp) also) 2=€ 4 y=. (7) 
eingefiihrt werden; es lautet dann: 
Ze, = GE, N, 2, 2, 2,) 


mit den Anfangswerten z(€, &), #(é, &), g(& &) auf der Anfangsgeraden 7=&. 


Wir schreiben im folgenden statt €, 1 wieder x, y, statt G wieder f, so daB 
also jetzt das Problem 


s=f(%,¥,2, 0,9), 
2(x,x)=a(x), p(x, *)=B(x),  9(%, x) =y(2) 


zur Lésung vorliegt. 


(3) 


Fiir die Entwicklung der numerischen Verfahren ist es niitzlich, die Gebiete 
y=x und yS<~ gesondert zu betrachten. Dabei soll die Gerade yx jeweils 
zu beiden Gebieten gehéren. Es werden daher ,,Formeln fiir Werte oberhalb 
y=x‘* und ,,Formeln fiir Werte unterhalb y= x‘ entwickelt. 

Sei (€, 7) ein Punkt oberhalb y= x, der zum Definitionsbereich des Problems 
(3) gehdrt, so ist dieses unter den Voraussetzungen des Satz 1 aquivalent den 
Integralformen: 


2(€,y) = 2, €) shal, x dx+f Sil 9,2, b.Q)aydx, 


ia esas +1 (é,y, 2.9) dy, (4) 


U] 


q(é,) =49 (n,n) a ee 2, p,q) ax. 
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Ist (€, 4) ein Punkt unterhalb y=, so wird entsprechend: 


2(n, 6) = 29,0) = Joe, Nat ee y,2,~,qQ) dy dx, 


n 


b(n, &) = p(y,.n) —StE%2%2,9) 4y, (5) 


4 
s 


a(n, &) = 9, €) + S10, z,p,q) 4x. 


Mit Hilfe der Streifenbedingung (2) laBt sich das Integral 


dia mee durch Kin fk dig lauds Pre = @iere (0) 
é é 


ersetzen und man erhalt gleichwertige Vorschriften. Es wird darauf aber nicht 
weiter eingegangen. 
Die Intervalle q<aS*<p<p,, y<ySyS6< 6, werden unterteilt durch 


Oe Neg ean 
ae (7) 
B= <M <0 <=. 
Es sei ferner stets 


Kuti Xp = Vv41— W=h, Fie ig Sek Bos Tee oma 8 


3. Extrapolationsverfahren 
Die Rechnung sei bereits so weit fortgeschritten, daB Naherungswerte 


ey 4+u,s—v> Pye Gr+pu,s—v Pea ea (7, S20; 4=0,1,...,m; v=0,1,. -., M— Uw), 


bekannt sind. Es sei bemerkt, daB die im folgenden verwendeten Approximationen 
nicht die einzig méglichen sind. Da durch Extrapolation aber hauptsdchlich 
Rohwerte gefunden werden sollen, diirfen wir uns beschranken. 

Es wird dann mit =x,, n=¥,4, in (4) f(x, y, 2, p, g) durch das Polynom 
ersetzt, das an den Stellen (%ypins Viae) = 04, ean nO et rea ledic 
Werte annimmt [1], [4]: 


IF, SS”. 
m m—pE . ) 
UM Se) LB *) 
2) 
mit 
hi hy SHAN, Gy VN Shs (9) 
Entsprechend wird ersetzt: 
m Ve» > 
f(%,, VrePs q) durch Sean Feiss 
v=0 : 
m Uu+ty—-1 
Bias Meaty Bo ‘De q) durch ar Ae beta st (10) 
u=0 
GL Gy 
q(x%,*+(y,—4%,)) durch > aI Ata, 
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Dabei ist gesetzt: 
U™) = u(w — 1) (uw — 2)...(u—w+1), 


VO" = ov +1) (v+2)...0+r—1), 


Tyw+1)—1 Ueé+)) 


U 


(Bei der Wahl =x,, n=y, 4, ist die Anfangsgerade y= x-+ (y,— %,).) 


Fiihrt man die Integrationen in (4) aus und beriicksichtigt, daB y= x-+ (y,— 


in v=u tbergeht, so erhalt man fiir Werte oberhalb y—x: 


m+1 mm m—U 
a 
ay s+1— *,5 ob 1)" Bi A’ g, — h? >) De pee Ved oes 
7] u=0 v=0 


Pr sti = rst hEP, Voi 15S? 


v=0 


m 


Gy, s+1 = Gy+1,s+1 alt Ds (= NE By Neiiey cia 
u=0 


lesa = f(%,, Vs+1> By s-+1> Pen: Upset) , 
mit 


Buv= 


U 


(—1) M+» 


uly! 
0 
1 
1 


fm 1)u(u+1)...(4+p—2) | f(0—1) 0-2)... (0-1) a du, 
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x,) 


(12) 


(13) 


6,= jefuoern. SG ANd uw, Bi=>, [ u—1)u(u+4)...(u+4—2) du 


0 


Tabelle 1. Gewrchte B, y: 


Bo,o =3% By,0 = sf Boo =—vF 
Boa =75 Bia = 22 
Bo, ee ve 


(Fiir die Gewichte 6, und fF siehe Tabellen 1 und 4 in Teil I.) 


Die Abschatzung der Kubatur- und Quadraturfehler wird bei allen Verfahren 
zur Lésung des Cauchy-Problems recht umstandlich. Da man bei der Extra- 
polation diesen Fehler wohl kaum beriicksichtigen wird, sei hier auf eine Ab- 


schatzung verzichtet. 


Es empfiehlt sich, die in (12) auftretenden Differenzen durch Funktionswerte 


zu ersetzen. Nach einiger Rechnung erhalt man: 


m+1 m m—lE 
es * 2 SA 
Sy s+] Sees ENO r a, 16 Opt den h Ds Loe CHEE en eee 
u=0 “u=0 v=0 


Peis baat bY tm sho vy? 


™ 


Gy, s+1 — Fr+1, POs whee tas s+1> 
pu=0 


if Sake = =} (%,, Vsti: Pperen s Py. st1> Yr, seas 
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mit a a 


Tabelle 2. Gewichte af ,: 


m=0: Ko =4 

m=1: et o= 8 at o=s 
o1= 8 

m= 2: to,9 =a ato = 39 oS o=— 37 
1=— Fe e1=—dz 
Oba = Fz 


(Zu den Gewichten a,,, und o,44 ,, siehe Tabellen 2 und § in Teil I.) 


Fiir die Entwicklung entsprechender Formeln fiir Werte unterhalb y= ~ sei 
vorausgesetzt, daB die Rechnung schon weiter fortgeschritten ist und Naherungs- 
Werle 2 ute n shies Prti—p,sti4or Ur ti—uctiew de tiaactd +e DOMOmnt sina (a 
1,...,m; v=0,1,...,m —). Es ergeben sich dann die (14) ganz ahnlich ge- 
bauten Extrapolationsformeln fiir Werte unterhalb y= x: 


m+1 m m—U 
a * SEG PRS NSORRES 
pe i= Ai ales ea i ee lO De ye eee 
u=0 u=0 v=0 


m 
Preto Pes iet h > Celi oer 
v=0 


(16) 


m™ 


Ah SS 
Gy +1,s cis Yr, s Ge h ves om u cae ) 
u=0 


h+1,s= carey Vs» Sy41,5> Preiss, FrAas) 


4. Berechnung des Anfangsfeldes 


Um die soeben in 3. und die in der nachsten Nummer 5. angegebenen Inter- 
polationsformeln fiir die fortlaufende Rechnung anwenden zu kénnen, bendtigt 
man schon entlang der Geraden y= vx ein ,,Anfangsfeld“ von Werten. Es soll 


hier speziell ein Interpolationsverfahren zur Berechnung dieses Anfangsfeldes 
angegeben werden. 


Auf der Geraden y=» sind die Anfangswerte z(x,, y,), D(%,. Vy), Q(%r> Vp) 
F(x,, y,) vorgegeben. Man verschafft sich dann zuerst Rohwerte z!),,, pl, ,, 
Gigs Polen, (7=0, 1 Se ee +m—r) mit 
Hilfe der Formeln (14) und (16), die etwas umgeschrieben werden und dann 
lauten: 
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Fiir Werte oberhalb y= x, (k=1): 


=i aly 
ae os = ye ith y oR Kine wre h 1—h?. > 2 Goes, y ns Peat y=1> 
i u=0 v=0 
oe = 3 ie eee +h 2 Ae—1» iy nae 
a (17) 
Cpe = dena, rp 23,1 (as 
ee 
es ee Vr4 my) Mk Be Lh» Gye a 
Fiir Werte unterhalb se a RSA) 
pal ay 
Ay = BY oe =e AD op dee eh a h? 4 a ack Fife ee sees 
be f r=0 
j= = pr} (ley ie her h x Xe—1 a rtk+l+y> 
aay (18) 
Gh = Oe gan th by ei 1 he i rie 
[0] 


[0 
[Shop aaa f(%,, Vrtko ® meee PrP Gay) . 


Dabei ist zu setzen: 


0 0 0 0 
a = 2% Ve)> Pr = P(X)» = Tr I)» =F (%,,)- 
Die so gewonnenen Rohwerte werden nun iterativ verbessert: 


setzt man einmal in (4) =%,, 7=y,4;, em andermal in (5) 7=4,, €=%y,45, 
so erhalt man — mit Beriicksichtigung der Streifenbedingung — nach ent- 
sprechender Approximation der Integranden durch eindeutige Polynome folgende 
Interpolationsformeln, die als Iterationsvorschriften benutzt werden: 


Bhren=2 (Xp Ie) thd ths m4 (Xu Vu) He] Saal (as Vu) mms Post » : Hes os is) 
u=0 


w= ar v=0 
M+ 


Beh P (Vy) +h | TP, mE (p> Yq) + > Me A Z| 
ae | (19) 


ae =% q(Xp+> a) Aa Hl ‘ F(x. ‘t+k? Vr +p) = Das ies 


1 1 
{Plea =f (Xs Veto ae srgs pe bigs qi bata) 


mit den Gewichten 


h 
—4)e+" : 
M2, = =e | git asain Har he 
0 
xX (utr—p—1)...(4utr—m) x 
7 (20) 
f@tn...@ Ly—y+1)(v+r—v—1)...(u+7r—m) dv| du, 
k 
k 
a Loe a (wty7)...(utr—A+1) (ut+r—A—1) ... (utr—m) du. 
! te 
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Tabelle 3. Gewichte nes 
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Die Formeln (19) gelten gleichermaBen fiir Werte ober- und unterhalb y=~x. 
Fiir A=1 vereinfachen sich die Vorschriften (19) zu 


Pees = ik a h® > Dons », mh (f14), — ised) ) 


u=0 v=0 
ptr 
p= Pe). + bY Way 
ee, (21) 
A+1 2 A A= 
qv jar ba — Daath Gea 
ce 


A+1 A+1 A+1 A+1 
en Ses Vrths ee py Ate ql oe 


Fir die exakte Lésung der Differentialgleichung gilt: 


z(%,, Vyik) = 2(x,, Vy) +h Sm m4 Ca Vedi 


“u=0 
+> Sf mF (%,, ie oe sh Bat ShIth, 
D(X. Vp ee) = p(%,, Vy) +h Sh F(a, 9) + SA 


Y (Xp. Vota) =U (%rte>» Veer) — Dt ee eval Ce Sosa 
Lt 
Fir die Kubatur- und Quadraturfehler seien hier nur einfache aber grébere 


Abschatzungen gegeben: 


gnrt ay ymrs k2 emrt emrl ui 
| m+, m+1| S m+1 2 || da”t1 we oym™tl ~ \Max 
mrt ozmr2 
' m+1 | ax™t1 Gyms wae 
= mre om 
[Sia] <| | sp me (23) 
m+1 | ex Max 
jmr2 om+1 
Si ale la — ye , 
m+1 | aymt Max 
mre om+l1 
Sit] <a as 
m+1 Bega Max 


Bei der Berechnung der Gewichte (20) ergeben sich Vereinfachungen: Wie 
man leicht bestatigt, ist [siehe Teil I dieser Arbeit S. 434]: 


(SOME OA (24) 
und fiir ¢=1 
in Nee (25) 
Ferner gilt nach (20): 
ee m == Aas = ak v,m ae ie m Ms Aek — hn) p (26) 


30* 
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SchlieBlich 1aBt sich die Berechnung der 7\”%,, noch vereinfachen: Ersetzt man 
in (20) nacheinander v durch —v+a, u durch —u+a, so wird 

a 
; | (u—a)(w—a+1)...d¢—a+tp—1)x 


ul y!(m—p)!(m—y)! 


xX (u—atw+ti)...(4a—atm) x (27) 


mE ae 


0,0 
Nu,v,m 


u 
x| f (v—a) (v—a+1)...(v—a+v—1) (v—a+yv+1)...(u—a+m)dv|du. 
0 
Mit Hilfe der Formeln (17), (18), (19), bzw. (24) errechnet man die Werte 
iiber dem in Abb. 1 (fiir m=2) mit ,,1“ bezeichneten Quadrat. Die gleichen 
ynr Formeln dienen aber auch zur Berechnung 
der Naherungswerte tiber den in Abb. 1 mit 
2‘, 3, ... bezeichneten Quadraten. Man 
braucht mit der Zahlung der Werte tber 
den einzelnen Quadraten lediglich wieder 
von vorn anzufangen. So wird man z.B. 
Sy ep anstelle VON Zig icc m-brag. SOvzeuy 
Cee ee 
Die Werte iiber den Quadraten stellen 
dann ein Anfangsfeld dar. Alle anderen 
Werte kann man schon mit dem Verfahren 
Yo x, Lp Ly Ls fay fiir die fortlaufende Rechnung ermitteln, 
Abb.1. Zur Berechnung des Anfangsfeldes allerdings unter Benutzung von héchstens 
10 Funktionswerten. Reicht die Genauigkeit 
nicht aus, so kann man sich mit der Berechnung weiterer Werte iiber ,,Zwischen- 
quadraten“ (in Abb. 1 schraffiert gezeichnet) helfen. Auch dazu kénnen offenbar 
die hier angegebenen Formeln dienen. 


5. Interpolationsformeln fiir die fortlaufende Rechnung 


Es werden wieder gesondert Formeln zur Berechnung von Naherungswerten 
ober- und unterhalb y= x entwickelt. Da die Betrachtungen hierfiir denen in 
3. sehr ahnlich verlaufen, kénnen wir uns kurz fassen. 


Mit £=*,, 7=y,,, wird in (4) ersetzt: 


SEM! Viet 
f(x, Y, 4, p, q) durch Dy b 7 i At Vi, le s+1> 
oe ob 
q(x,x+(y,—-%,)) durch Sat An 
u=0 - 
™m (W141)? (28) 
1 (G5 2sP, @) durch s * Vedra 
»=0 
m ul) 
f(*, Vsi19%, p, q) durch pu! dn i 
=0 ; 
Dabei ist ; 
UM = u(u —1)...(¥u—p+1), 
(29) 


(V— 41)" = (© = 1)o@+ 4)... (toa), 
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Fiihrt man dann in (4) die Integrationen aus und ersetzt die Differenzen durch 
Funktionswerte, so erhalt man als Iterationsvorschriften die Interpolations- 
formeln: 


m m—pE 

[A+1] __ 2 am {V4 

By, Sail aay arth a dete, Sty 1 Xo, Eerie, er 1 RyS4+ lv | 9 
pu=0 u=0 »v=0 


1, v= 0,0 


m 
A+1]__ * A * 
Meow: (Dee + h Chea ae “hs See coe , 
v=1 


(30) 
aif * fA * 
Tee a Ble =k Loge © rt+p,st1|> 
oa 
nie) Al), AMG) Ieren 
Med = 1 Op, Vente Oy te Geel) 
mit 
mM—v M—HK x” y 
aot ae ap * 
Spiel Ae oy Dt)e eal lear 
== “AP, 
* (ee ; 
£.=4SG fe t)ee +1)... @+m—2) x (31) 
0 
i 
x f 2@—1)...@—9+1)de 
0 
ae , wie (15). 
Tabelle 4. Gewichte Bz ,: 
Bo,0=% Bio=s 2,0= —1'r 
Bo.1=—% Bi1=—y¥z 
Bo, 2= — 12 
Tabelle 5. Gewichte aj"): 
m=0: eo o=F 
mat: ah=t ahh 
og, 1=6 
M=2: 09,985 ot o= 28 3 = — 34 
9 12 of ive 
09,2= —¥E 


Fiir die Berechnung der Naherungen unterhalb y= x erhalt man entsprechend 
die Gleichungen: 

m s m mM" a 
==/}) pe: Coe UPR aS eee pti a Ds Cary Tot Lee ASate Pal ee 


u=0 u=0 y=0 
bv 0,0 


ress 
Ei sey ay oad 


D 


m 
2 * 
= Py 4 1,s+1 —h lado fet > Am, v lease 
r= (32) 


m 
naa cag Dah ahr , 


A A+1 
pire (Spas Vs» gly pt, a) 3 


A+1 
Mi=q th 
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Fiir 421 vereinfachen sich die Iterationsvorschriften betrachtlich: 


Aus (30) wird: 


A pg k= = BE ae a Wao. 6° 0 ( ale = ea) . 
Se = Pea re od ee 
Ge oh =O a Ny oe at Nae lg 
He = fps Segip oe ty Peel) eo): 
Aus (32) wird: 

ay ees = By, s — hase 0 ( ee Feat Ma; 
= pl, hone ie — HD), 

qt = gis, oa han, olen, s Sire © 
aE HCO y Senso eben eerie 


(34) 


Fir die exakte Losung der Differentialgleichung gelten ebenfalls die Gleichungen 
(30) und (32), wenn dort anstelle der Naherungswerte die exakten Werte eingesetzt 


werden, mit den Kubatur- und 


Quadraturfehlern (s. (22)) aS en 


m+1 ae 1St3 


See ise, fiir Werte oberhalb y= x und *S7 44 mia + aS as 25% LL 25% | 1 fiir 
Werte unterhalb y=. Es gelten die Abschatzungen: 


| oy eae 


m+1 


| s < pmt3 ba om. 
2 m+1—A95,,A ore 
| Shen m+1 y= Ox ay Max 
15% 
| m+1 em 
2 = Soy ee Zea : =H Beer: 
Sed Max 
Cx 
| Spee a ymre omri ‘|B: 
Cc: = aym+1- m+1 
OF oe, | Oy Ma 
1S* 
m+1 om+1 
m-+2 oO * 
ra Sh are earns 
25% Ox Max 
m+1 


Die Fehler der Naherungen fiir z, die nur durch genaherte Kubatur und 
Quadratur entstehen, haben die Form 


Pek, wee es oe 
Ss A hm? as By — Chpmt?. 
fess stato <a lessee 


Formal erreicht man Proportionalitat der Fehler zu h”*3, wenn 


m+1 

in (30) Ss pe A es qd, + ob, Se 
ss 

in (32) 


u=0 


m 


statt > Opn, 0 Vy +u, Soe? 


u=0 


> hi yu Vp l—pg,s+1—y statt S che lu VE ajo ein bey 


u=0 


gesetzt wird. 


Es kénnen noch andere Interpolationsverfahren gefunden werden, je nach der 


Art der Approximation der Integranden in (4) und (5). Prinzipiell ergeben sich 
aber wohl keine neuen Gesichtspunkte. 
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Es sei schlieBlich noch erwahnt: 

Kennt man bei der Behandlung des Cauchy-Problems schon eine Anzahl von 
Naherungswerten an bestimmten Stellen, so kann mit den Verfahren zur Lésung 
des charakteristischen Anfangswertproblems weitergerechnet werden (siehe Teil I). 


6. Konvergenz der Anfangsiteration 


Fiir die folgenden Betrachtungen geniige das Problem (1) den Voraussetzungen 
des Satz 1. Setzt man 


a a 
[opel + ope + | oy3| = ef, 
so gelten fiir die ae ae, (21) die Abschatzungen 
of Bove aM | 3 Ingle & inch wil, 
u=0 »v=0 v=0 u=0 
+r Ver wmer+k (37) 
C=O at RSS, HT. oc, — 1, 1. me 7 
Dieses Ungleichungssystem hat die Form 
Ce Sa, > ee elke a, (38) 
=O) y=) 
; Mev 
Sind alle Eigenwerte der Matrix ., - 
B= (Bi) 
dem Betrag nach <1, so ist die Konvergenz der Reihen 
P Oo”! y+k» y > 64 rt+k» oy OV. (39) 
=o = 
gesichert. Hinreichende Bedingungen fiir die Konvergenz sind 
Max 13 > Sait <1,  (Zeilensummenkriterium) , 
a2=0 750 
Hore (40) 
Max » yi oe. lr ,  (Spaltensummenkriterium). 
Hwy? |y=0 k=—r 
k+0 
Sei m 
Max > Sule ad ? Max >). |7¢2,(=— By Max > es =, 
te y= 0 v=0 1% y=0 me u=0 
bey vaeY uer+k 
so ist nach (37), (38) und (40) sicher 
hW27MA+hM(B+C)<1 (41) 


hinreichend fiir die Konvergenz. Sei ferner 


Max {A, B, C}=L 
dann ist Se 
A 
n<|At ap eo 


ein weiteres hinreichendes Konvergenzkriterium. 
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Fiir m=2 lauten die Bedingungen: 


Nach (40): Rae Neonat 
<8 irr) 
aOR 
Seale Y 375M ay 
Nach (41): ts Bick (43) 
25 ny 
oe pee +r 4). 


Nach (42): 


7. Konvergenz der Iterationen bei der fortlaufenden Rechnung 
Mit den Bezeichnungen (36) gelten fiir die Gleichungen (33) die Abschatzungen : 


ws (12M |ag o| + 24M |o%, ol) oH. (44) 


Ein hinreichendes Konvergenzkriterium fiir das Verfahren (33) bzw. (30) und 
— wie man sofort iibersieht — auch fiir (32) bzw. (34) ist daher: 


M (i? |ak | 2h +|oX ol) <1. (45) 
Hieraus erhalt man noch 
h< cee oe ae =ih (46) 
0! Tl (Cee 0) y 
Fir m=2 ist 
n<to(|/1+58 —1). (47) 


8. Fehlerabschatzung fiir die Naherungswerte des Anfangsfeldes 


Es werde vorausgesetzt, daB die Iterationen (21) jeweils bis zum Stillstand 
durchgefiihrt wurden, so daB die Indizes [A] fortgelassen werden kénnen. Ab- 
rundungsfehler werden nicht beriicksichtigt. 


Der Einfachheit halber sei der Quadraturfehler Stertk in in (22) zum Kubatur- 


fehler Si,"4¥41 hinzugenommen, so daB im folgenden statt S%"V",,44-+ Saara nur 


sie nya Zeschrieben wird. Subtrahiert man die Gleichungen (22) von den Glei- 
chungen (19), so gelten mit den Bezeichnungen 


Buy cz Z (Xu yy) = Eu Bis a p (%» yy) aa ae Guy > q (Xu, yy) iF EM, (48) 


die Abschatzungen : 


lenenl FAM DDI ml {leat [ell €u,o|}= [Seas 
: Mey 

eee —hM la mile Ey} |S) y| + eo tt ss | Sak , (49) 
Ve 

le, re ~ Hat Il Pie Chern eh oat 


ae 
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mit den Abschatzungen (23) fiir die Kubatur- und Quadraturfehler. Die Matrix 
des Systems (49) ist von monotoner Art [5], mit 


|e, Nes Za |é, “alae |é, p= |Z as (50) 
Site als SSEIBL ST yc Combed Re Bae ripe oy iva bese Pte 
gilt daher: 
(Zs 1k eM 2 2 > \nes v, bal y+Z,, eZ, MUS i bein 
a 
Zope — bMS (nhl {Zt Ziv + Z}= Belt, (51) 
er 
Le a hu > Gee WNZE: r+k +Z, tk aL I, sean Te 
see 
Das Gleichungssystem (51) hat die Form 
je oe 


und mit der Einheitsmatrix © die nach (40) konvergente Matrizenreihe 


s= (CF BLVL--)F 


als Lésung. Sei 


Mas ie lta mae Shll —a. 


Asay. 
so gilt bei geniigend kleinem / fiir das gréBte Element Max (8) der Matrix 8: 


Max (8) =hMa 
und allgemein fiir das groBte Element Max ($”) der Matrix %”: 
Max (8”) = h” M"a A*-1, 


Fiir die Naherungswerte des Anfangsfeldes gilt dann die Fehlerabschatzung: 


k Mha Tu, = i 
lé,r+e1S Sears prey oe sy AC mt, mde Ly Ty) , 


—hMA 
207 = 
bv 
Panis Ra ay hMa aa Tie Te Ti 2 
| Seg | = se Te ye ( ed, wa wet ae a (5 ) 
u=0 v=0 
MY 
Trrtk hMa a. Te To? Tio 
| Sea m+1 Lear eeg eV aT: DURE Bia eli 
pu=0 v=0 
we 


Damit ist die Fehlerabschatzung im Prinzip geleistet. 
Mit Hilfe von (52) laBt sich zeigen: Zu beliebig vorgegebenem e>0 1aBt sich 
stets eine Schrittweite # so angeben, daB fiir 7=0, 1, ...,.m; k=—yr7, —r+1, 


1 1 U7 


Nees i os teem, |e; eee ete: 
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9. Fehlerbetrachtungen beim Interpolationsverfahren der fortlaufenden 
Rechnung 
Es seien die Voraussetzungen aus 8. erfiillt und obere Schranken der Fehler- 
betrige des Anfangsfeldes bekannt. Die folgenden Fehlerbetrachtungen werden 
fiir die Formeln (30) zur Berechnung von Naherungswerten oberhalb y= x durch- 
gefiihrt. Fiir unterhalb y= x gilt Entsprechendes. 
Mit den Bezeichnungen (48) gelten dann fiir die Fehler die Abschatzungen: 


le, saa SL igils Joa’ | {le 54a 5 eee 1| =a Koen Auge 
a hM |op,ol{|& sxa| 37 (ep elles [&, s+} S Am» (53) 
oP Es” bi 


|é 


erat 
7,s+1 


ae 1| =F AM |omn,ol{|é,s+al + ees nai 


mit 


Ay,=|é r, [+43, coma oes oye 


tee x 
u=0 v=0 
4,v+=0,0 


x1 Sus tes a | Seip teaser trate [eren aimee te | tS ats l= ceed) 


By = lé,, sick ia lomo {| 2, eae ae ES etes| sh [Pe Sa eyieee eco , 


ark, =| Gacy cal bh 3 | om, ul {le 1+ th, Peat\ te | eae seal tléey, ealt+hs 15|- 


Die A,,, Ap; A,, sind nach Voraussetzung bekannt, ihre Berechnung allerdings ist 
miihsam. Es lassen sich aber obere Schranken angeben, die einfacher berechnet 
werden kénnen: 


Sei 


Max {| See eee Maxtlen | rasa yeoz | Max {|@,44,041-/}=Z, 
(U0, 4, 24.905 We 0,45. ptt nee = 0, Os 


ferner 
m mL m 
dX Llane |=Gim)=G, DY |oe,,|=L(m) =L, 
u=0 vp=0 wal 
4, v= 0,0 


so gelten die Abschatzungen 


Ape Ey, Ata (L+ lok |Z) 4+WM G{Z4+Z4Z}+|} Ss ile She ese 
Ay, S|, | HAML{Z+Z4Z}+| 5% |=Ry, (55) 
Ae ele e: s41| FAM L{Z4+Z4Z}+ 1S* SI=R,, . 


Unter Beriicksichtigung der Monotonie der Matrix des Systems (53) erhalt man 
dann die Fehlerabschatzungen: 


h? M |ok% Dailies 
le ace aw ie, Re a7) —h? M |\at™| — 2h M[ox, al (Ry +R, + R,)> 


[Estill [Rn h Max, o| R +R S6) 
Eval) [Ry SP MTeoR ahaa, | Bm 7 


ee 
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Damit ist gezeigt, daB wenigstens prinzipiell eine Fehlerabschaitzung moglich ist, 
wenn obere Schranken der Fehlerbetrage vorher berechneter Naherungen bekannt 
sind. 


10. Ein Konvergenzsatz 
Es wird im folgenden gezeigt, daB bei der numerischen Behandlung des Cauchy- 
Problems die mit dem Interpolationsverfahren fiir die fortlaufende Rechnung 
gewonnenen Naherungswerte bei kleiner werdender Schrittweite gegen die Werte 
der exakten Lésung des Problems (3) konvergieren. Fiir die Extrapolations- 
verfahren laBt sich der Beweis ganz ahnlich fiihren. 


Es sei das Anfangsfeld berechnet, (eventuell mit einem anderen Verfahren 
als dem in dieser Arbeit beschriebenem), d.h. es sind Naherungswerte 


ate, Atrth> Phtr,aAteth> Gate, Aretk: (57) 
(A=0, m, 2m, ...; 7=0,1,...,05 R=—7, —7 +1, ..., =1, +1,...,m—7), 
bekannt. Die zugehérigen Fehlerbetrage seien 
leraea4 pte |> Fear, Pater 9 [ete ae | (58) 


ferner 


Max{|é,4,,247+2|} =Zo> Max {|é,4,,.+r+0/}=2o> Mand | een iran} Lo: (59) 
Es werde nun vorausgesetzt : 


lim E, = lim £, = lim E,=0. (60) 


h—0 h—>0 h—>0 


SchlieBlich seien die in den Abschatzungen fiir die Kubatur- und Quadraturfehler 
vorkommenden Ableitungen hdéherer Ordnung von F(x, y) beschrankt. Dann 
gilt der 


Satz 2. Ber der numerischen Behandlung des Cauchy-Problems seven 
By s> prs Vp, 5 (2 Bok aes S == 0, 55 anos 7=<S),, 


die mit dem Interpolationsverfahren der fortlaufenden Rechnung ermittelten Nahe- 
rungswerte der exakten Losungswerte 
2(%p, Vs), Ply Ws)» 1% Vs) 


des Problems (3). Dann ist bei beliebigen aber festen a, 6, wenn z(a, B), p(x, B), 
q(«, 8) die exakten Werte und 2,3, bx, p> Iu,p de Naherungswerte sind: 


2a, g — (a, B) = &,,5> 0, Pup — P (a, B) = €x,6 +0, Gap — 7 (¢ B) = 8,2 > 0 
fir h—>O. 


Beweis. Es werden alle Naherungen der z, f, g, die langs der Geraden y= 
xt (u+1)h, (u=1, 2, 3,...), mit dem Interpolationsverfahren fiir die fortlaufende 
Rechnung ermittelt werden — die also nicht zum Anfangsfeld gehoren — mit 
Zu» Pur %, bezeichnet. Die zugehérigen Fehlerbetrage seien eA sleniene 


E,|- 
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Sei Ey le, 
Xo ar Ey , bu as Ey, ? 
E é 
ferner : : (61) 
IS. 
34 = gk: | Sul 
|S. 


der Vektor, dessen Komponenten die maximalen Betrage aller Kubatur- und 
Quadraturfehler sind, die bei der Berechnung der z,, ,, 7, langs der Geraden 
y=x+(u+1)h entstehen, so gilt, wie man leicht erkennt, fiir das Interpolations- 
verfahren der fortlaufenden Rechnung eine Abschatzung der Form: 


W 41S Biko +8, =t, (62) 
mit 
‘By (h)  By(h) by 3(h) 
B= C+) bs). beelh) dosh 1. 0-0, ab fe ba, Via 
bsi (A) b32(h) 033 (h) 
ha, ha, ha, 
Me OR as de Gaal MO Rajge ies eee Deu: (64) 
ay ag Ag 


Die Matrix Y% ist monoton, denn es existieren xy, und der vom Nullvektor ver- 
schiedene Vektor r>0, so daB 
Sei 


so folgt wegen der Monotonie von Y aus (63) und (65) fiir h=0: 
ace tye calsOy ether l ya hake \é, |< Ey. 
Nach (65) ist 


Ue = Ato + (BV — A) ty +4. (66) 
Da alle Elemente der Matrix § — Yl fiir h>0 positiv sind, so folgt: 
AQ@j92ae, unddaher ea, 
wegen der Monotonie von XA. 


Das soeben beschriebene Verfahren zur Ermittlung oberer Fehlerschranken 


kann fortgesetzt werden. Allgemein gilt fiir die Fehler der Naherungen langs 
der Geraden y= «+ (n+1)h: 


Wen = Bt + 3y, (=AV 2,805), (67) 
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mit 


v 


Ln = cin Opi Op. aes lg/Sz, 


~ 


| &} & 


n 


Wird &4B=C gesetzt, so folgt aus (67), wie man durch vollstandige Induktion 
bestatigt: 


n—1 
GO tor DO U4 sg, mit C=C. (68) 
pu=0 
Die Vektoren %18,,_,, miissen Darstellungen der Form 
0, —,(h) 
ae pen —— Meas Oyen (h) 
6, —,(h) 


besitzen mit 


ee) eS eee enrs (25,95 2 1G ya Ar}| < BOM eattlTes ren 
Die Matrix % hat nach (64) die Form 
Y= €—Y. 
Die Matrizenreihe 


konvergiert, wenn simtliche Eigenwerte der Matrix &% dem Betrag nach <1 
sind, was nach (45) stets erfiillt ist. Sie stellt dann die Matrix 


yA (E- A) 
dar. Die Matrix € kann dargestellt werden durch 
Cy1(2) Cro (A) C13 (h) 
C= CAP Hint) Cy5(2) “Colwet,  O<¢,.< C. 
C31(2) C3a(h) C33 (h) 
Bezeichnet man das gréBte Element der Matrix ©” mit Max (©), so gilt: 
Max (8) < 1+4+AC 
Max (©?) < (4+4C)(1+3h0C), 


Max (6) <(1-+hC) (14340 < (14 3AC)4. 


Betrachtet man jetzt die Fehlermajoranten auf den festen Geraden y= x + k mit 
(n+1)h=k, so gilt mit nh=k—hA<K fiir alle wSn: 


Max (6) < (14+3A4C)%"<e8°* fir ho. (69) 
Nach (68) folgt daher 


by 


n 


E,\ <&*(E, +E, +EothttK(S+5+S) (70) 
E, 
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mit 
Max{S,_,}=S, Max{S,_,}=S, Max{S,_,}=S, (w=1,2,....%—1). 
Wegen der Beschranktheit der héheren Ableitungen von F(x, y) sind auch die 


S, S, S beschrankt. Nach den Voraussetzungen am Anfang dieser Nummer folgt 
dann sogleich ae. = 
lin Eee 0 jodi BOs la ea) 
h—>0 h—>0 h—>0 
und daher auch p ges ges 
lim |2,(2=0,- lim \eghe= 0, sim [elie 
h—>0 h—>0 a) 
und daraus die Behauptung des Satzes. 
Bei der Extrapolation ist &%=%-!= € und somit C=%. Der Beweis 1aBt sich 


hier fast wortlich genau so fithren wie fiir die Interpolation. 


11. Beispiel 
Differentialgleichung : 


Org, Ge. 02, A PY 
ax Oy Ox By | fag Fs wan a 2 to 
Anfangsbedingungen: 


2 (ae) =e" cosas 
Dix, ie 3 4.6." COS, 
q(x, x) = e?* (xcos x — sin x). 
Losung: 
ge Kegs ys 
Ergebnis der Rechnung: 

Das Anfangsfeld (umrandet) wurde mit m=2, h=0,1 nach (17), (18) und 
(19) bzw. (21) berechnet, die fortlaufende Rechnung mit Verfahren (30) bzw. 
(33) und (32) bzw. (34) mit m=2, h=0, 1 durchgefiihrt. Bei der Anfangsiteration 
und bei der fortlaufenden Rechnung waren zwei bis vier Iterationen erforderlich. 
(Die Schrittweite / ist hier schon reichlich groB bemessen.) 

An den Randern des Gebietes wurden folgende Formeln benutzt: 

1. Oberhalb y= x: 

[eee 1 
By, 8 = Fy gy — 420 (91,6 = V,2— 59r43,8 + 0,05 13 at 0,35f,,9 72 
~~ OATS. Si 0,35 f,41,8 aia 0,05 f,41,7 sad 0,1 f,+2, 8) 


eh 1 
20.3 = CN aires (6907 + 6413+ 0,05 fas + 0,35 fo.7— 0,4 fo,6+ 


2. Untethalb y= x: + 0,35 f1,3 + 0,05 f1,2 — 0,1 fo,8) - 

44,0 = 843.1 — oo (59,0 + 8p41,1 — Frna,a + 0,05 (0 + 035 faa 1 — 
CLR pata = 0,4 fy aie 10:35 Goce O,05:f, 054g eee) 
ethan a (692.1 + 697, + 0,05 /8,b + 0,35 fs,1— 0,1 fe 2 + 


+ 0,35 17,6 + 9,05 f2.1 — 04 fo 0)- 
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